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ABSTRACT 

Chaff is the least expensive and simplest soft kill measure 
against anti-ship missiles. In this thesis we study and simulate 
the deployment of chaff in centroid mode for anti-ship missile 
defense. The idea is to fire chaff of radar cross section greater 
than that of the ship, cause a shift of centroid and achieve break 
lock by moving the ship out of the radar resolution cell of the 
missile. After firing the chaff in an optimum direction, the ship 
is made to maneuver in two steps; first to make an earliest 
possible break lock and then to maximize its displacement from the 
missile line of flight. Since the interaction scenario is time 
variant, modelling has been structured on a relative coordinate 
system. The algorithms to implement these objectives are based on 
values extracted from this coordinate system centered on the chaff, 
mobile with the wind and oriented towards the missile. 

The mathematical model is translated into a simulation 
program. We consider it to be a good tool for insight and 
understanding and feel that it has a promise of application both 


1n design work and at sea in the passive point defense problem. 
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THESIS DISCLAIMER 
The reader is cautioned that computer programs developed in 
this research may not have been exercised for all cases of 
interest. While every effort has been made, within the time 
available, to ensure that the programs are free of computational 
and logic errors, they cannot be considered validated. Any 
application of these programs without additional verification is 


at the risk of the user. 
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I. INTRODUCTION 


A. IMPORTANCE OF DEFENSE AGAINST ANTI-SHIP MISSILES 

The maritime combat scenario of today is filled with complex 
multi-dimensional threats for a man of war. Conflicts in the past 
quarter century or so have shown that warships have grown more 
vulnerable to aerial threats, in particular to anti-ship missiles. 
We have the following prominent episodes of recent times to endorse 
and highlight this fact: 


e An Israeli destroyer EILAT was sunk near Port Said in the 1967 
Arab-Israeli war by Egypt with Styx missile fired by Soviet built 
Komar class boats. 

e The Pakistan Navy destroyer KHAIBAR and minesweeper MUHAFIZ 
were sunk and logistic ship DACCA was damaged by salvoes of Styx 
missiles fired by Ossa class missile boats, which are also Soviet 
built, during the 1971 Indo-Pakistan war. 

e Gabriel missiles launched by the Israeli Navy inflicted heavy 
damage to the Syrian Navy resulting in sinking of a number of 
missile boats during the 1973 Arab-Israeli war. 

e The Royal Navy suffered a humiliating blow when a number of 
its front line ships were sunk and damaged by shore and air 
launched Exocet missiles fired by the Argentinian Navy during the 
Falklands war in 1982 (Ref. 10). 

e US ship and air launched Harpoon missiles destroyed a number 
of Libyan missile boats and other craft in March 1986 during a 
conflict with Libya in the Gulf of Sidra. 

e In April 1988 again, Harpoon and Standard missiles fired by 
the US Navy caused heavy damage to Iranian Navy vessels in the 
Persian Gulf. 


B- APPROACHES TO SHIP DEFENSE 

All major and minor maritime forces of the world recognize the 
vulnerability of surface ships to this serious threat today. 
However, since the surface fleets are foreseen to be the mainstay 


of most navies for a long time to come, their defense against anti- 


ship missiles has received tremendous emphasis. Larger navies who 


have access to sophisticated technological and vast monetary 
resources tie their belief to the methodology of layered defense 
or defense in depth. Such a defense methodology would probably 
include offensive long-range aerial strikes against enemy storage 
depots and platforms, combat air patrols, standoff missiles and a 
variety of other active or passive tactics. Even these elaborate 
measures fail to provide an absolute defense. Smaller navies who 
do not have such luxuries available at their disposal must live 
with lesser and cheaper methods. They have to rely on early 
detection of the enemy and in most cases early detection of the 
missile. Such a restriction imposes the use of short-range active 
or passive measures. These measures include active ECM, 
conventional guns, point defense systems or close in weapon systems 
(CIWS) such as Vulcan Phalanx/Breda and chaff. Ideally one could 
achieve most effective results by having all of them working in 
unison in an intelligent manner. However, availability of active 
ECM and CIWS involves sizable finances and in many cases their 
availability is dependent upon politics. Chaff, on the other hand, 


presents itself to be a cheap and simple alternative. 


Ce MODES OF CHAFF DEPLOYMENT 

Chaff has been used effectively in a number of military 
conflicts since the second World War. Its flexibility enables it 
to be used at various stages of conflict in many modes. The 


discussion here will be restricted to a confrontation between an 


anti-ship missile and a surface unit at sea. Chaff, in such a 
situation can be used in the following modes (Ref. 13:p. 338]: 
ho. Confusion 
In confusion mode, chaff is fired at long range (a few 
kms) when the ship is outside the range of enemy search and target 
indication radars. This tactic provides a large number of false 
targets that are indistinguishable from the real ship. 
2. Distraction/Dilution 
When the ship has been detected by the enemy and a 
missile attack is imminent, chaff is fired at shorter ranges (1-2 
kms) in order to mislead the missile radar when it starts 
transmitting. The chaff cloud's radar cross section need only be 
comparable to the ship in this mode. 
3. Dump 
This mode is also called range gate pull off/stealing. 
Once the missile is locked on, chaff is fired at short ranges (less 
than 1 km). In this mode the ship's ECM equipment is also used in 
conjunction with the chaff. The ship's jammer misleads the missile 
on to the chaff and when jamming ceases, the missile homes onto the 
chaff decoy. In other words, break lock is achieved with the help 
of the jammer. 
4. Centroid 
In this mode, chaff is deployed at very close ranges. 
The exact figure depends on the philosophy of actual deployment. 
Some experts suggest firing chaff at a range of 25 to 50 m so that 


its echo merges with that of the ship and the missile homes onto 


the missile is lured away avoiding a direct hit thus minimizing 
damage. This is a last ditch method as break lock is not the 


primary goal. 


D. RECOMMENDED SOLUTION 

Although reaction times are shrinking with sophisticated 
missiles entering the battle arena, intelligent deployment of chaff 
has been known to present an effective counter. Ideally, chaff 
could be most effective if integrated with other forms of active 
and passive anti-ship missile defenses. However, if monetary 
constraints preclude such a luxury, we consider that chaff provides 
a good alternative to thwart most missile threats of today. It can 
be argued that recent advances in technology have enabled 
production of smart and sophisticated missile radars which have the 
ability to reject chaff clutter by sensing the difference of its 
parameters from those of the ship, thereby rendering it 
ineffective. Nevertheless, there is always another counter to a 
counter. Chaff systems are also getting smarter with added 
features. These features include chaff’s maneuvering capability 
with the help of a self propulsion system inside the chaff and 
expendable ECM devices with integrated processing facility. Active 
decoy developments such as Nulka and Winnin are examples of such 


chaffs. 


Es AIM OF THE THESIS AND RELATED WORK 
The purpose of this work is to simulate centroid mode of chaff 


deployment, fire the chaff at close range (100 to 400 m), achieve 


break lock and maximize miss distance by intelligent ship maneuvers 
assuming a supporting jammer is not available. Work has also been 
done in this field by David Nitzan and Chris Nygard [Refs. 22 & 23) 
The mathematical modelling and simulation done in this thesis has 


been developed as an extension to their work. 


II. KINEMATICS OF THE PROBLEM 


Aa MISSILE DETECTION 

One of the most formidable and decisive problems in the field 
of anti-missile defense is early detection of the missile. Most of 
today’s anti-ship missiles travel at speeds close to MACH 1 and fly 
very low to the surface of the sea. This sea skimming capability 
coupled with the state-of-the-art aerodynamic characteristics, 
makes the task of missile detection even more difficult for the 
Operations Room of the ship. For vessels operating without any air 
support, reliance has to be placed on sensors available on board. 
The two key sensors of help in this situation are radar and ESM 
equipment. ESM can be coupled with radar to confirm detection made 
difficult by background clutter at low altitudes. However, modern 
missiles commence radiating for their terminal homing phase quite 
late in their flight, thus limiting the ESM capability to ranges 
less than the radiation commencement zone. 

The primary sensor, as a matter of principle, should be radar. 
Nonetheless, because of low flying profiles of the missile and non- 
avallability of state-of-the-art radars to smaller navies, ESM 
takes over the primary sensory role. This is especially true for 
the sea skimmers. The range at which the missile commences its 
terminal phase has a very large variability, depending on the type 
and generation of the missile technology. Based on available 
unclassified data, this range varies from around 3 miles to 20 


miles. With an average speed of 600 knots, at which most missiles 


travel, we get a reaction time ranging from 18 seconds to about 2 
minutes. It is quite obvious that reaction time is the key factor 
in today's warfare. With the missiles becoming supersonic and their 
terminal phases rapidly shortening, reaction time is increasingly 
diminishing. Table 1 shows reaction times for different initial 
detection ranges and is a useful indicator to emphasize this point. 
Assuming the missile speed to be 600 knots, the following data has 
been tabulated: 
TABLE 1 


MISSILE RANGE AND REACTION TIME 


Y 
zZ 
O 


Missile Range (yds) Reaction Time (secs) 


20,000 
15,000 
14,000 
13,000 
12,000 
11.000 
10,000 
9000 
8000 
5000 
3000 
1000 TOO LATE 


l. 
2. 
Br, 
4. 
at 
6. 
7. 
8. 
9. 





B. CHAFF DEPLOVMENT IN CENTROID MODE 

There are two philosophies associated with centroid or 
seduction mode tactic; the first is to fire the chaff at very close 
ranges, i.e., 25 to 50 meters (Ref. 5], within the radar resolution 


cell of the missile. The missile takes the two echoes to be one and 


locks on to the centroid of the two. The ship then moves away 
briskly from the chaff, depending on the wind direction leaving the 
missile to home on the chaff. Depending on all the variables 
involved, the ship can either break lock if it is out of the 
resolution cell in time or else can maximize displacement from the 
missile impact point. In this way, the ship can either escape 
damage altogether or at least, minimize it. Real time data from 
relevant conflicts is not available to seek a recommended figure 
for safe distance of the ship from the point of impact. However, 
at least 50 to 100 meters is considered to be a safe distance for 
minimal damage to the ship. 

The second philosophy behind utilization of chaff in centroid 
mode is to fire the chaff within the resolution cell on a bearing 
related to the wind direction and speed. Chaff is placed as far as 
possible from the ship towards the boundaries of the resolution 
cell [Ref. 13), depending on the initial detection range of the 
missile, beam width of the missile radar and the wind parameters. 
Once the chaff has bloomed, the missile amends its homing on 
direction to the centroid. In the meantime, concurrent to chaff 
launch, the ship maneuvers its way out of the resolution cell to 
achieve a break lock and maximize distance from the missile impact 


point. Analysis of this philosophy is the aim of this thesis. 


Es RADAR RESOLUTION CELL 
The resolution cell is the central factor of the anti-ship 


missile problem. Basically, the radar has a certain range and 


bearing resolution based on its effective pulse width and beam 
width (Ref. 3:p. 332). The area outside these bounds is the limit 
of the radar to distinguish one target from another. If there are 
two targets and the radar is locked onto one of them, the second 
should be outside the boundary to be discriminated as a separate 
target. If the second target also lies within this area, the radar 
will only see them as one target and align itself to their 
centroid, both radially and azimuthally. Figure 1 shows what a 


radar resolution cell looks like. 


Range Gate 
Beamwidth O +) te 





Missile 





ACTUAL 
Angle Gate 


Figure 1. Radar Resolution Cell [Ref 23] 

The limits in both angle and range are defined by the half 
power or 3 dB values. Most missile seekers use this geometry for 
processing purposes. From the missile's point of view, only the RF 
energy inside this cell is considered to be the useful backscatter. 
Any return from the target outside this area is ignored to ensure 
sharp resolution. The cell actually assumes an annular shape, but 


at ranges we are interested in, it can be approximated to be a 


rectangle without incurring any sizable inaccuracies as shown in 


Figure 2: 


Beamwidth 


Missile 
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Figure 2. Approximation to the Resolution Cell (Ref. 23] 
The bounds of the cell are also called gates; the one 
perpendicular to missile direction is the angle gate and the one 


parallel to it is the range gate as shown in Figure 3. 


MISSILE 


RANGE GATE 


q 
ANGLE GATE 


Figure 3. Range and Angle Gate 
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The center of the cell is the centroid for the missile radar 
to lock on. As the bearing of the missile changes due to movement 
of the ship or centroid, the cell rotates to reorient itself in 
order to maintain orthogonality with the missile line of flight. 
Similarly, due to variation of target range from the missile as it 
closes in, the dimensions of the cell vary in the azimuthal plane. 
The radar resolution cell, therefore, forms an orthogonal 
coordinate system which is varying with time and constantly 
readjusting its bounds. This motion, as mentioned, depends on the 


movement of the centroid and the longitudinal axis of the missile. 


D. EFFECT OF CHAFF ON RADAR RESOLUTION CELL 

The whole idea of using chaff in a favorable manner is to 
shift the centroid of the resolution cell away from the ship, walk 
it off to the chaff and move out of the cell to break lock thereby 
seducing the missile onto the chaff cloud. The chaff, of course, 
has to be fired to bloom inside the resolution cell. If the radar 
cross section of the chaff is comparable to that of the ship, the 
missile is expected to shift its lock to the centroid of the two 
radar cross sections which is approximately the mid-point on the 
line joining the ship and the chaff [Ref. 13:p. 347). 

Usually, most captains would like to fire chaff with a radar 
cross section considerably larger than that presented by the ship. 
The shift of the centroid would then be more favorable to the ship 
and closer to the chaff cloud. In any case, as soon as the chaff 


blooms, it starts moving with the wind. The aim is to fire the 
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chaff on a bearing that optimally assists the ship for the wind 
vector existing during interaction. The job of the ship then is to 
move out of the resolution cell, in a direction that will take it 
away from the chaff. 

The objective of moving the ship out of the cell is to reduce 
and then eliminate its radar cross section from the bounds of the 
cell. As the ship starts moving away from the chaff, its distance 
from the centroid also increases. As it crosses the cell boundary, 
either through its range or angle gate, its radar cross section 
inside the cell decreases. The effect of this development is two- 
fold: first the ship has moved away from the centroid; second, due 
to its decreased radar cross section inside the cell, the centroid 
moves further towards the chaff. This effect multiplies as more of 
the ship exits the cell and hastens the shift of centroid to the 
chaff. As the ship leaves the cell completely, break lock is 
achieved. The missile now locks onto the chaff and the ship 
endeavours to maximize its distance from the new centroid, i.e., 
the chaff. 

The kinematics are, however, not as simple as they appear to 
be. The most crucial factor at this stage is the wind speed and 
direction in relation to both the missile line of flight and the 
ship's initial course. The wind could be a helping hand or prove 
unfriendly depending on what it is doing during the operation. If 
it is strong and moving parallel to the missile line of flight, 
i.e., perpendicular to angle gate, the chaff will spend very little 


time in the cell. This time might not be enough to enable the ship 


12 


to exit the cell before the chaff does so. As a result, the missile 
locks on to the ship again. 

As an exception, the ship and the chaff might decide to leave 
the boundaries of the cell simultaneously. In such a situation, it 
is hard to tell what the missile will lock on to. There appears to 
be an equal probability of the missile either locking on to the 
ship or the decoy. However, it first depends on the decision-making 
algorithm of the processor of the missile's radar. It would also 
depend on the relative location of the chaff and the ship with 
respect to the missile. If the processor unit conducts a far-to- 
near processing pattern, the nearer of the two will be selected. 
Of course, the reverse would be true for a near-to-far pattern. 

In summary, the presence of chaff redefines the boundaries of 
the radar resolution cell, reorients it, displaces it and aids the 
ship to get out of troubled waters. The crucial factor is a 


favorable wind vector at the scene of the battle. 


je 


III. MATHEMATICAL MODELLING 


In order to formulate a simulation program, it is necessary 
to construct a mathematical model of the situation in which the 
chaff is to be deployed. There are a number of parameters that are 
made available to us by the ship, missile, wind and the chaff. Some 
of these parameters are constant and others are variable. The 
variables are classified into two categories; those that are 
presented to us by a particular situation or the input variables, 
and those that we have the ability to control and manipulate. It 
is the latter type of variables that make it possible for us to 
successfully use the chaff. A list of input and control variables 
are given in Tables 2 and 3. The abbreviations for all the 
variables in this chapter are chosen to correspond to those used 


in the simulation program. 


A. SHIP*S EQUATIONS 

As summarized in the tables, the ship's parameters consist of 
its length, breadth, course, speed, acceleration, turn rate and 
radar cross section. 

je Initial Conditions 

The ship is assumed to be at the center of the reference 

coordinate system for the initial set up at the time the missile 
is detected. This is actually the screen center. The initial 


velocity vector of the ship is determined by its initial course and 
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speed.The following equations describe the ship 


point: 
XS (h) 
YS (h) 
XSD(h) 


YSD(h) 


XCENT 
YCENT 
VS (h) *Sin{ (CSE(h) *7/180} 


VS (h) *Cos( (CSE (h) *7/180) 


TABLE 2 


LIST OF INPUT VARIABLES 


VARIABLE 


Initial Bearing 

Initial Range 

Missile Speed 

Radar Beamwidth 

Range Gate 

Missile Acceleration 
Proportional Navigation Constant 


Course 

Speed 

Length 

Breadth 

Turn Rate 
Acceleration(Velocity Rate) 
Max Radar Cross Section 
Min Radar Cross Section 


Radius 
Radar Cross Section 


Wind Direction 
Wind Speed 


MISCELLANEOUS 


20. 


Screen Scale 


r5 


at the starting 


(seal) 
(2702) 
(3.3) 


(3.4) 


ABBREVIATION 


BRGM (1) 
RANGE 
VMM 

6 
RGATE 
AM 

N 


CSE (1) 


RCSSMAX 
RCSSMIN 


CRAD 
RESE 


GW 
VW 





TABLE 3 


CONTROL VARIABLES 


VARIABLE ABBREVIATION 


First Leg of New Course CSE2 
Second Leg of New Course CSE3 
Max Speed VSMAX 


Bearing to launch chaff BRGC 





where 

XS(h) and YS(h) are coordinates of ship displacement 

XSD(h) and YSD(h) are coordinates of ship velocity 

XCENT and YCENT are coordinates of screen center 

CSE(h) is the ship's course 

2 Ship’s Movement 

Having detected the missile, the ship after firing the 

chaff, accelerates and maneuvers out of the radar resolution cell. 
In order to update its position, velocity and acceleration 


components, the following equations are used: 


XS(h+1) = XS(h) + XSD(h) + 0.5 XSDD(h) *h? (3.5) 

YS(h+1) = YS(h) + YSD(h) + 0.5 YSDD(h) *h? (3.6) 
where 

XSD(h) = VS*Sin{CSE(h)*r/180) (3:44) 

YSD(h) = VS*Cos(CSE(h)*r/180) (3.8) 


The velocity components of the ship are shown in Figure 4. 
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e CSE(h) VS(1) 


XSD(h) 


Figure 4. Ship's Velocity. Components 


Also 
VS(h+1) = VS(h) + VR*h (3.9) 
XSDD(h) = -VS(h)(TR*x/180) * U * Sin(CSE(h)*x/180) 
+ VR*Sin(CSE(h)*1/180) (3.10) 
YSDD(h) = -VS(h)(TR*r/180) * U * Cos(CSE(h)*r/180) 


+ VR(h)*Cos(TR*r/180) (3 11) 
We have supposed that VS can attain a maximum value of 25 
knots. The terms used in the above equations are defined as 
follows: 
XSDD(h) and YSDD(h) = Instantaneous coordinates of ship's 
acceleration 


XSD(h) and YSD(h) = Instantaneous coordinates of ship's 


velocity 

VS = Velocity of ship 

VR = Velocity rate or acceleration of the ship 
U = Turn direction (discussed later) 

h = Time step increment 
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Updating the bearing of the missile from the ship is carried 


out using the following equations: 
BRGM(h) = 180/r * tan” E htl} —_¥s (nt) | (3510) 


where XM(h+1) and YM(h+1) are updated coordinates of the missile 
position. 
3: Radar Cross Section 
The radar cross section of the ship is a variable 
dependent upon the aspect presented by the ship to the missile 
radar. It is assumed that the missile presents maximum radar cross 
section when viewed from right ahead, right astern and the beams. 
As the aspect changes, the radar cross section decreases to a 
minimum at the bows and the quarters. In the absence of real time 
data, we consider this assumption to be a reasonable approximation. 
It is modelled with the help of the following equation: 
RCSS = RCSSMIN+ |RCSSMAX-RCSSMIN*Cos[2(CSE(1)-BRGM(1))*A/180]| 
(3.13) 
where 
RCSS = Radar cross section of the ship 
RCSSMAX = Maximum radar cross section of the ship 
RCSSMIN = Minimum radar cross section of the ship 
CSE(1) = Initial course of the ship 


BRGM(1) = Initial bearing of the missile 
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Pictorially, the equation produces a graph of the ship's radar 


cross section vs the aspect angle as shown in Figure 5. 


000 


270 090 


225 135 


180 


Figure 5. Radar Cross Section Model 


ES 


4. The Ship's Model 


The ship's geometry is approximated by a polygon centered 
on (XS,YS) using its length and breadth. The orientation of the 
ship, pivoted on its center is shown with the help of its course 
with respect to true North as shown in the Figure 6. 

The equations to find its vertices are given as under: 


XSHIP(1) = XS(h) + SL/2 * Sin (CSE(h)) 


(3.14) 
YSHIP(1) = YS(h) + SL/2 * Cos (CSE(h)) (3.15) 
SHIPS MODEL 
North 
CSE 





Xship(5), Yship(5) Xship(1), Yship(1) 


Xship(2), Yship(2) 


Figure 6. The Ship's Model [Ref. 23] 
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XSHIP(2) = XS(h) + SL/4 * Sin (CSE(h)) + SB/2 * Cos(CSE(h)) 


(3.16) 





YSHIP(2) = YS(h) + SL/4 * Cos (CSE(h)) - SB/2 * Sin (CSE(h)) 
| (3.17) 
XSHIP(3) = XS(h) - SL/2 * Sin (CSE(h)) + SB/2 * Cos (CSE(h)) 
(3.18) 

YSHIP(3) = YS(h) - SL/2 * Cos (CSE(h)) - SB/2 * Sin (CSE(h)) 

(3.19) 
XSHIP(4) = XS(h) - SL/2 * Sin (CSE(h)) - SB/2 * Cos (CSE(h)) 
(3.20) 

YSHIP(4) = YS(h) - SL/2 * Cos (CSE(h)) + SB/2 * Sin (CSE(h)) 
(3:21) 

XSHIP(5) = XS(h) + SL/4 * Sin (CSE(h)) - SB/2 * Cos (CSE(h)) 
(3.22) 

YSHIP(5) = YS(h) + SL/4 * Cos (CSE(h)) + SB/2 * Sin (CSE(h)) 


(3:28) 


B: CHAFF EQUATIONS 
The following parameters are related to the chaff: 


Launch bearing 

Distance from the ship 

Radius of the chaff 

Radar cross section of the chaff 
Bloom time 

Velocitv vector of chaff after bloom 


E Initial Conditions 
Since the wind is the fundamental issue here, it is 
appropriate to see its relationship with the chaff. Wind direction 


is, by convention, taken to be the direction it is blowing from, 


Sch 


the direction in which the chaff moves is its reciprocal. This is 
shown in Figure 7. The unscaled velocity of the wind is VW1 and it 


is also the velocity of the chaff. 


Velocity Vector 
North 


Wind 






GC=GW+180 


VW1=VC 


Figure 7. The Wind Velocity Vector 
where 
VW1 = Wind speed (Unscaled) 
GW = Wind direction 


Velocity vector of the chaff is given by the following 


equation: 
VC = VW 
GC = GW + 180 (3.24) 


Let us now consider the initial positioning of the chaff. In 
order to offset the centroid as far away as possible, we fire the 
chaff so that its center lies on the boundary of the angle gate. 
The angle gate is, however, shrinking as the missile closes in. 


The initial placement of the chaff, therefore, is brought closer 
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to the cell center laterally by one chaff radius in order to make 
it bloom inside the future dimensions of the cell. The calculation 
of this distance is done using the following equation: 

CDIST = R(1) * TAN(THETA * PI/180) - CRAD (3.25) 


The initial position of the chaff round is shown in Figure 8. 


Launching Vector 


North 


CDIST Cos BRGC 





Ship 
wi CDIST Sin BRGC 
(XCENT, YCENT) = (XS, YS) 


Figure 8. Initial Position of Chaff [Ref. 23] 
The parameters in the figure are as follows: 


BRGC 


Chaff launch bearing 
CDIST = Distance of chaff from the ship 
We can see here that the coordinates of the chaff are: 


XC(h) = XCENT + CDIST * Sin(BRGC * 2/180) (3.26) 


YC(h) = YCENT + CDIST * Cos(BRGC * 1/180) (3.27) 
2. Chaff Movement 
The ship is initially in the screen center. As the chaff 


starts blowing with the wind, its coordinates are updated by the 


following equations: 
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XC(h+1) = XC(h) + VC * Sin (GC * r/180) (3.28) 


YC(h+1) 


tl 


YC(h) + VC * Cos (GC * 1/180) (3.29) 
3. Bloom Time Delay 

In order to Maximize the factor of realism, we have taken 
into account the time delay between the firing of the chaff and its 
full bloom. For simplicity, it has been assumed that the time taken 
to launch the chaff is included in the bloom time delay. Although 
various types of chaff have different time delays in this context, 
we consider a bloom time delay of 5 seconds as an appropriate 
average figure (Ref. 13:p. 342). It also includes the time it takes 
to launch the chaff at its intended initial position. The 
implmentation of this concept is achieved by keeping the chaff on 
the ship for the first 5 iterations. On the sixth second, the chaff 
cloud comes into existence at the desired location. Thereafter it 
is at the mercy of the wind vector and starts drifting. Since the 
ship and the resolution cell are mobile, there is a possibility 
that the chaff may bloom outside the cell as a result of bloom time 
delay. To avoid such an eventuality, the first position of the 
bloomed chaff is offset in accordance with the ship's movement. 

4. The Chaff Model 

The chaff cloud is approximated by a circle but due to 
convenience of mathematical calculations, it is drawn as a hexagon 
using graphics. The hexagon is centered on the chaff center (XC,YC) 
and is shown in Figure 9. Coordinates for the corners are worked 


out with the help of the following equations: 
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XCF (1) 


YCF (1) 





Figure 9. Chaff Cloud 
XCF(2) = XC(h) + CRAD * Cos(60º 
YCF(2) = YC(h) + CRAD * Sin(60° 
XCF(3) = XC(h) - CRAD * Cos(60º 
YCF(3) = YC(h) + CRAD * Sin(60' 
XCF(4) = XC(h) - CRAD 
YCF(4) = YC(h) 
XCF(5) = XC(h) - CRAD * Cos(60º 
YCF(5) = YC(h) - CRAD * Sin(60° 
XCF(6) = XC(h) + CRAD * Cos(60º 
YCF(6) = YC(h) - CRAD * Sin(60° 

©. EQUIVALENT TARGET EQUATIONS 


I 


XC(h) + CRAD 


YC(h) 


XCF(3), YCF(3) 


XCF(5), YCF(5) 


XCF(2), YCF(2) 


XCF(1), YCF(1) 


XCF(6), YCF(6) 
Model [Ref. 23] 


) 


Se tn aap 


w 


(3.30) 


(3.31) 
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Before the chaff is fired, the missile radar is Locke on to 


the ship and is tracking it. When the chaff is fired and it has 


bloomed inside the radar resolution cell of the missile, the radar 
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senses a new situation. The weighted center of the two, i.e., the 
ship and the chaff, is now seen by the missile as the centroid of 
the target and it starts tracking it. This centroid, also called 
the equivalent target, is determined by the magnitude of the radar 
cross sections of the chaff and the ship as seen by the missile. 


It is calculated by the following equations: 


XET(h+1) = XS (h) *RCSS(h = XC(h) * RCSC(h (3.42) 
YET(h+1) = YS (h i— h - YC(h) * RCSC(h (3.43) 


The velocity with which this equivalent target moves is calculated 


using the following equations: 


VET 


[(XET(h) - XET(h-1))“ + (YET(h) - YET(h-1)}°]® (3.44) 


= = YET(DN) = El 


This is shown in Figure 10. 


D. MISSILE EQUATIONS 
The missile parameters consist of the following: 


Range and bearing from the ship 
Range resolution 

Missile speed 

Half power beamwidth of the antenna 
Proportional navigation constant 
Acceleration of the missile 

Turning rate 
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(XET(h), YET(h)) 


{XET(h-1)}, YET(h-1)} 


Figure 10. Equivalent Target Velocity 
l. Initial Conditions 
The first parameters related to the missile that are 
observed by the ship are its initial detection range and bearing 
which are shown in Figure 11. The related equations are as follows: 
XM(i) = XCENT + R(i)* Sin(BRGM(i) * 7/180) (3.46) 


YCENT + R(i)* Cos(BRGM(i) * 1/180) (3.47) 


YM(i) 
The direction flight of the missile is the reciprocal of its 
bearing from the ship i.e., 

GM = BRGM(i) + 180 (3.48) 
The initial velocity components of the missile in the X and Y 


directions are shown in Figure 12. The equations for these 


components are as follows: 


XMD(h) = VM * Sin(GM * 2/180) (3.49) 


il 


YMD(h) = VM * Cos(GM * 1/180) (3.50) 
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North 
Missile 


(XM, YM) 
BRGM(h) 


R(h) Cos BRGM(h) 
(XCENT, YCENT) 
R(h) Sin BRGM(h) 


Figure 11. Missile Range and Bearing 


North 


XMD(h) GM(h) 


(XM, YM) 
YMD(h) 


Figure 12. Missile Velocity Components 
25 Missile Movement 
The coordinates of the missile position as it moves are 


updated as follows: 


XM(h+1) = XM(h) + XMD(h) * h + 0.5 * XMDD(h) * h° (3.51) 
YM(h+1) = YM(h) + YMD(h) * h + 0.5 * YMDD(h) * h° (3.52) 
XMD and YMD are extracted from the equations in the previous 
section whereas, XMDD and YMDD, which represent the components of 
the missile acceleration are shown in Figure 13 and given by the 


the following equations: 


XMDD = AM * Sin(GM) (3.53) 
YMDD = AM * Cos(GM) (3.54) 
North 
Missile 
(XM, YM) 
GM 
AM YMDD 


(XCENT, YCENT) 
XMDD 


Figure 13: Missile Acceleration Components 
The missile direction angle is updated using the following 


equation: 
= sl XM(h+1) - XM(h 
p- | M(hr1) — YM(N) | E-P) 
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This is shown in Figure 14. The X and Y components of the missile 


position vector are updated iteratively. 


(Xm(h+1), Ym(h+1)) 







YM(h+1) - YM(h) 


(XM(h), YM(h)) 


XM(h+1) - XM(h) 


Figure 14. Updating Missile Direction 
SS Target Tracking by the Missile 
In order to maintain track of the target the missile uses 
the principle of proportional navigation constant i.e it maintains 
a constant line of sight with respect to a reference. This angle 
is shown in Figure 15, 


where 


= d YET(h)_- YM(h 
o = tan | ERM | (3.56) 


In order to maintain a constant line of sight angle, the 
missile guidance system senses the rate of change in o. This is 


given by the following equation: 
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RMT (Ch) 
(YET-YM) 


REFERENCE 
(XET-XM) 


Figure 15. Target-Missile Line of Sight (Ref. 23] 
where 
RMT(h) = Range of the missile from equivalent target 
In order to cater for this rate of change in the line of 
sight and bring it to zero, the missile needs to turn at a 
corresponding rate. This is called the turn rate or the 
acceleration of the missile and is given by the following equation: 
AM = SIGD * VM * n (3.58) 
where 


n = Missile navigation constant 


E.: RELATIVE COORDINATE SYSTEM 
T. Shift of Origin 
In the complex and randomly occurring scenario of 
different ship headings, missile bearings and wind directions, 
movement of the chaff after full bloom can take many paths in terms 


of relative motion. In order to minimize the disadvantageous 
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features of the variables involved, we have decided to have a 
relative coordinate system centered on the chaff round. This is 


accomplished with the help of the following equations: 


XR(h) = XS(h) - XC(h) (3.59) 


YR(h) 


YS(h) - YC(h) (3.60) 
where XR(N) and YR(N) are the relative abscissae and ordinates 
respectively. A sketch depicting the shift of the reference system 


is shown in Figure 16. 


SHIP 
All 
(A) 
CHAFF (NEW ORIGIN) 
(2.1) 
(0,0) 
OLD ORIGIN 
(0,0) 
(= 2 out ) 


Figure 16. Shift of Reference System 
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The original coordinates of the positions of the ship and chaff 
are (1,2) and (2,1) respectively but when we shift the origin to 
the chaff, ship's position using the above mentioned equations 
becomes (-1,1). 
2. Rotation of the System 

The next step is to rotate the new coordinate system to 
the bearing of the missile in the clockwise direction. Such a 
rotation is desirable and handy to calculate the chaff launch 
bearings which are orthogonal to missile direction. This is done 


using the following relationship: 


| Cos BRGM(h) -Sin BRGM(h) 
XRP(h) = | | XR(h) 
| Sin BRGM(h) Cos BRGM(h) | (3.61) 


where XRP and XR are the new and previous state vectors 
respectively. Upon multiplication with the matrix we get the 


following equations: 


XRP(h) = XR(h) * Cos BRGM(h) - YR(h) * Sin BRGM(h) (3.62) 


YRP(h) 


ti 


XR(h) * Sin BRGM(h) + YR(h) * Cos BRGM(h) (3.63) 
This mathematical model can be better understood using Figure 17. 
We see that point (0,2) becomes (-1,1.7) when the system is 
rotated. 
oie Relative Velocity Components 

The next important step is to calculate the relative 
velocity components using the initial values of the ship’s velocity 
and the velocity of the wind with which the chaff moves. This is 


shown in Figure 18. 
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Figure 17. Rotated Reference System 


Using Figure 18 we get the following equations: 


XRPD(h) = VS*Sin(CSE(h)-BRGM(h)) + VW*Sin(GW - BRGM(h)) (3.64) 

YRPD(h) = VS*Cos(CSE(h)-BRGM(N)) + VW*Cos(GW - BRGM(h)) (3.65) 
where VS and VW are the velocities of ship and wind respectively 
and XRDP(h)/YRDP(h) are the X and Y components of the relative 


velocity respectively. These values of relative velocity components 


are calculated using the initial conditions. 
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Figure 18. Relative Velocity Components 


Is. PRE-COMPUTATION METHOD 

The random variabilitv of the missile direction, ship's course 
and speed at the time of missile detection along with wind speed 
and direction, present a much more complex problem than is readilv 
obvious. It is extremelv difficult to obtain the best solution in 
all the possible cases with conventional chaff alone. Efforts are, 
however, made in this work to enhance optimization. In order to 


achieve this, it is considered useful to carrv out a rough pre- 
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computation which can give an idea as to what chaff firing bearing 
would entail the largest miss distance. Making use of the relative 
coordinate system, the problem becomes somewhat simpler in 
formulation. As mentioned earlier, the chaff is fired orthogonal 
to the missile bearing, leaving us with two alternatives;one to the 
right or in the positive XRP direction and the other in the 
negative XRP direction i.e., to the left. The distance traveled by 
the chaff in either direction is computed for the entire run of the 
play. A couple of assumptions made in this regard make this a 
rather rough but a pretty good estimate of the situation. The first 
assumption is that the ship's speed is taken to be 25 knots which 
is the maximum speed of the ship assumed in this work rather than 
a gradual increase from the ship's initial speed. The next 
assumption ignores the alteration to the first leg of the new 
courses which is carried out to escape the range gate in minimal 
time. A direct alteration to the second leg from the ship's initial 
course is assumed instead, which is perpendicular to the missile 
direction. The following equations calculate the magnitude of the 
relative motion of the ship with respect to the chaff cloud in 


either direction. 


XRP1 = XRP1 + XRPD1*T + 0.5*XRPDD*T* (3.66) 

YRP1 = YRP1 + YRPD1*T + 0.5*YRPDD*T* (3.67) 

XRP2 = -XRP1 + XRPD1*T + 0.5*XRPDD*T* (3.68) 

YRP2 = -YRP1 + YRPD1*T + 0.5*YRPDD*T* (3.69) 
where 

XRPD1 = XRPD1 + XRPDD*T (3570) 
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YRPD1 = YRPD1 + YRPDD*T (3:71) 
and 


XRPDD 


-ACC*Sin(CSEDIF1) (3.72) 


YRPDD = -ACC*Cos(CSEDIF1) (3.73) 
where 
XRP1/YRP1 = Distance traveled by chaff in X/Y direction 
XRPD1/YRPD1 = Velocity component of chaff in X/Y direction 
XRPDD/YRPDD = Acceleration component of chaff in X/Y direction 
ACC = Acceleration of the ship 
CSEDIF1 = Difference in the ship“s initial course and the new 


course 


This is shown in Figure 19. 


G. RANGES AND BEARINGS 
1. Ranges 
There are three different ranges of interest shown in 
Figure 20 that are required for subseguent calculations. The range 
between the missile and the ship center is given by: 

R(h) = [ (XS(h) - XM(h))“ + (YS(h) - YM(h))* 15 (3.73) 
Similarly the range between the missile and the chaff center is 
given by: 

RCM(h) = [ (XC(h) - XM(h) )? + (YC(h) - YM(h))* J* (3.74) 
and the range between the missile and the centroid or the 
equivalent target is: 


RMT(h) = [ (XET(h) - XM(h))% + (YET(h) - YM(h))* )* (3.75) 
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Figure 19. Pre-Computed Distances 
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Figure 20. Various Ranges from the Missile 
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2 Bearings 


Like the ranges mentioned in the previous section, there 
are also some bearings that are necessary for subseguent 
calculations. These bearings are shown in Figure 21. The bearing 


of the eguivalent target from the missile is given by: 


BRGMT(h) = tan“ | | (3.76) 


NORTH (XS,YS) 


BRGMS (XC,YC) 
CHAFF 


RS, 
S 


(XM,YM) 


Figure 21. Various Bearings from the Missile 


3.9 


The reciprocal of this bearing, i.e., the bearing of the missile 


from the equivalent target is 


BRGTM(h) = tan“ | esp HE | (3.77) 


The bearing of the ship from the missile is: 


BRGMS (h) = tan’! | dr, ré | (3.78) 


Lastly, the bearing of the chaff from the missile is: 


BRCMS (h) = tan”! | ir re | (3.79) 


H. RANGE AND ANGLE GATES 

The range gate of the missile is approximately a function of 
its pulse width and is an input variable to the problem. The angle 
gate, on the other hand, is dependent upon the distance of the 
missile from the target and the half power beamwidth of the missile 
radar antenna. This is calculated using the following equation: 

AG = RMT(h) * Tan(6,) (3.80) 
We also need to calculate the distances of the ship and the chaff 
from the center of the radar resolution cell, which is the 
equivalent target or the centroid of the two radar cross sections. 
The position of the ship with respect to the radar resolution cell 
center along the range and angle gate directions, is calculated as 


follows: 
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AGS = R(h) * Sin | BRGMS(h) - BRGMT(h) | (3:81) 

RGS = | RMT(h) - R(h) * Cos(BRGMS(h) - BRGMT(h)) | (3.82) 
This is shown in Figure 22. Similarly, the position of the chaff 
is given by the following eguations: 


AGC = RMC(h) * Sin | BRGMC(h) - BRGMT(h) | (3.83) 
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Missile 
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Figure 22. Ship“s Position Relative to Cell Center 
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RGC = | RMT(h) - RMC(h) * Cos(BRGMC(h) - BRGMT(h)) | (3.84) 


This is shown in Figure 23. 


I. RADAR RESOLUTION CELL MODEL 
The radar resolution cell is constructed with the help of the 


range and angle gates of the missile radar. It is centered on the 


Chaff 
(XC, YC) 
AGC 
Resolution 
Cell 
RGC N 
RMC Equivalent Target 
A YET) 
a RMT 
GMT 
(XM, YM) 
Missile 


Figure 23. Chaff Position Relative to Cell Center 
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equivalent or the centroid and is approximated to be a rectangle 
as discussed in Chapter II. The cell is oriented towards the 
missile bearing. As mentioned earlier, the range gate of the 
missile is constant. The angle gate, on the other hand, is 
dependent on the half power beamwidth of the antenna and range of 
the missile so it varies as the missile closes in. The cell is, 
therefore, a time variant rectangle increasing in size as the 
missile“s range decreases and reorienting to maintain orthogonality 
with the missile direction. 

The depth of the resolution cell is twice that of the range 
gate and the width is twice that of the angle gate. This is shown 


in Figure 24. 
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Figure 24. Radar Resolution Cell Model (Ref. 23) 
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The edges along the width of the center line are calculated 


as follows: 


XW1 = XET(h) + AG * Cos((BRGM(h)) (3.85) 
YW1 = YET(h) - AG * Sin{ (BRGM(h) } (3.86) 
XW2 = XET(h) - AG * Cos((BRGM(h)) (3.87) 
YW2 = YET(h) + AG * Sin((BRGM(h)) (3.88) 


The corners of the cell are calculated as follows: 


XCELL(1) = XW1 + RG * Sin((BRGM(h)) (3.89) 
YCELL(1) = YW1 + RG * Cos((BRGM(h)) (3.90) 
XCELL(2) = XW1 - RG * Sin( (BRGM(h)) (3.91) 
YCELL(2) = YW1 - RG * Cos( (BRGM(h)) (3.92) 
XCELL(3) = XW2 - RG * Sin((BRGM(h)) (3.93) 
YCELL(3) = YW2 - RG * Cos((BRGM(h)) (3.94) 
XCELL(4) = XW2 + RG * Sin((BRGM(h)) (3.95) 
YCELL(1) = YW2 + RG * Cos((BRGM(h)) (3.96) 
J. PORTION OF THE SHIP INSIDE THE RESOLUTION CELL 


Je Ship Fully Inside the Cell 

The entire radar cross section of the ship lies inside 
the resolution cell initially. After firing the chaff, the ship 
attempts to maneuver out of the cell. In the process, the center 
of the cell shifts to the equivalent target. The new situation 
alters the radar cross section of the ship as seem by the missile 
in order to calculate the centroid. This cross section depends on 
the portion of the ship inside the cell and is a function of the 


following parameters: 
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Radar resolution cell center 

Size of the cell 

Relative location of the ship center 
Ship's heading 


The RCS of the ship is assumed to be the maximum presented by 
the ship at a particular aspect angle when it is fully within the 


bounds of the cell as shown in Figure 25. 


Resolution Cell 


te 





Figure 25. Limit for Maximum RCS of the Ship 
in the Resolution Cell (Ref. 23) 


2 Ship Partially Inside the Cell 
When some part of the ship is outside the cell, we have 
a partial radar cross section case. In order to calculate it, we 
need to shift our reference system from true North to the 
horizontal. This shift is necessitated by the fact that the 
computer calculates angles from the horizontal reference as shown 


in Figure 26. 
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Figure 26. Angle Translation to Fortran [Ref. 23] 


A conversion is, therefore, undertaken to cater for this new 


reference system. We define a new angle B, as shown in Figure 27 
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Figure 27. True North to Fortran Conversion (Ref. 23] 
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between the new reference and the direction from which the missile 
is approaching the ship. The aspect angle of the ship is called 
“Saa” and it is calculated using B. It is defined as the angle 


between the ship's heading and B as shown in Figure 28. 





Resolution 
Ship Center Cell 
XS, YS 
SB 
Range 
Gate 


Beta 


Figure 28. Aspect Angle [Ref. 23] 

The ship iS approximated as a rectangle for ease of 
calculation. Also shown in this figure are two more parameters; 
namely Slr and S2r. The subscript ’r’ merely represents the fact 
that the ship is exiting through the range gate boundary of the 
cell. If it exits via the angle gate boundary, then we replace it 
with subscript ’a’. S2r is the shortest distance between the range 
gate boundary and the ship center and it is given by: 

S2r = | RG - RGS | (3.97) 
Slr represents the portion of half of the ship’s length that is 
inside the cell and is calculated with the help of the following 


equation: 
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S2r 


Sir = os (Saa 


(3.98) 


In the case shown in Figure 28, where the ship center is 


inside the cell, the portion of the ship's RCS in the cell is given 


by: 
PARTS = — + 82-51 (3.99) 
where 


S1 = Sir or Sla depending upon how the ship exits the cell 
Sarea = SL*SB 
In case the ship center lies outside the resolution cell, the 


calculation is as follows: 





PARTS = —} - 22:31 (3.100) 


If the ship is exiting a corner (leaving the range and angle 
gate boundaries simultaneously), then the amount of ship's RCS in 
the cell is approximated by using the lesser of the two values 


obtained from Sir or Sla. 


K- PORTION OF CHAFF INSIDE THE CELL 
1. Chaff Fully Inside the Cell 

The chaff is fired to bloom inside the radar resolution 

cell. After it blooms the wind takes over and the chaff cloud 


starts moving with the wind vector. Since the cell is centered on 


48 


the equivalent target, some portion of the chaff might, dependent 


on 


the situation, spill out of the cell. In order to get an 


estimate of the location of the center of the cell, it is important 


to 


know how much of chaff is inside the cell. This depends on the 


following parameters: 


of 


of 


Relative location of chaff center 
Size of chaff 
Location of resolution cell center 
Size of the cell 
The chaff lies completelv within the cell when the distance 
its center from the cell bounds is at least equal to the radius 


the chaff. This is shown in Figure 29. 
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Figure 29. Limit of Maximum Chaff RCS in the Cell (Ref. 23] 
2. Chaff Partially Inside the Cell: 


We define the following parameters to assist in 


calculation of the chaff’s RCS inside the cell: 


Carea = 7 * (CRAD)“ (3.1010 
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where 
CRAD = Chaff radius 


| RGC - RG | (3:102) 


Rir 


Ria | AGC - AG | (3.103) 
If the chaff is located at the range gate boundary, the 
situation will be as shown in Figure 30. In this particular case, 


more than half the chaff lies outside the bounds of the cell. 





Figure 30. Portion of Chatt RCS in the Cell (Ref. 23] 
If we expand the view of interest, as shown in Figure 31, we can 
calculate a few more parameters that would assist us in finding 
chaff RCS inside the cell. 


We see that: 





= Si Rir 
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The sector bounded by twice this angle is given by: 


SECTr = 1/2 (2*A2r) CRAD* (106) 








Chaff within the 
Resolution Cell 


Gate radius 


Chaff Cloud 


Figure 31. Partial Chaff RCS in the Cell Blown Up [Ref. 23] 
Now 
R2r = CRAD * Sin(A2r) (3.106) 
The area of the two triangles, which is also the portion of the 
sector outside the radar resolution cell is: 
TRIr = 2(1/2*Rlr*R2r) = Rlr $ R2r (3107) 
In case the chaff exits the cell through the angle gate, 
the same procedure is applied using subscript Rla. If both the 
gates are cut off simultaneously, then we use both procedures to 
calculate the area or the chaff RCS inside the cell. All three 
possible cases are shown in Figure 32. 
When Rir or Ria are zero individually, the portion inside 


the cell is one half of the total. However, when they are both zero 


bi 


simultaneously, as in the corner case, then it is only one quarter 


of the total. 
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Figure 32. Various Cases of Partial Chaff RCS [Ref. 23] 
If RG >RGC, then most of the chaff cloud lies inside the 
resolution cell. The portion of chaff in the cell is calculated as 


follows: 


2 CAREA - SECTr + TRIr 


IF RG < RGC, then the RCS in the cell is less than half and is 


given by: 
CS ECTS RER ET 


When the chaff exits the cell from a corner we have four 
situations: 
1. RGC < RG and AGC <AG 
This is shown in the Figure 33 along with the actual and 


approximated radar cross sections. 
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Figure 33. Chaff Exiting the Cell Corner (Ref. 23] 


RGC < RG and AGC > AG 


2. 


This situation is shown 


Figure 34. 
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3. RGC> RG and AGC < AG. This is shown in Figure 35. 
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Figure 35. RGC.GT.RG £ AGC.LT.AG [Ref. 23] 
4. RGC > RG and AGC > AG 
In this case the larger of the two areas produced in the cell by 
RGC or AGC is approximated to be the chaff RCS of interest as shown 
in Figure 36. Less than a quarter of the chaff is inside the cell 
and it is taken to be one half the area cut off by the lesser AGC 


or RGC. 
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Figure 36. RGC.GT.RG € AGC.GT.AG (Ref. 23] 
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L. CHAFF LAUNCH BEARING AND SHIP MANEUVERS 
L. General Description 
In order to maximize time to break lock and miss 
distance, a double maneuver by the target ship is simulated. The 
first maneuver is to get out of the radar resolution cell of the 
missile via the shortest route, i.e, either towards the missile or 
directly away from it. The next maneuver is to alter to a course 
perpendicular to the missile direction when break lock is achieved 
and the ship is outside the resolution cell of the missile’s radar. 
2. Chaff Launch Bearing 
First of all, the bearing to launch the chaff is worked 
out. It is basically desired to launch the chaff at right angles 
to the missile direction. The two possible true bearings are 
calculated as follows: 


BRGC1 


BRGM(1) + 90 (3.110) 


BRGC2 


BRGM(1) - 90 (3.111) 

The decision regarding the chaff firing bearing is 
dependent upon the results of the pre-computation in order to 
optimize the miss distance. If the ship travels more along the 
positive XRP direction, we fire the chaff towards BRGC1 and vice 
versa. 

In the particular case shown in Figure 37, BRGC is equal 
to BRGC2. It must be kept in mind that the missile is shown to come 
from the top because of the rotated reference system and could in 


actuality be approaching from any direction in the azimuth. 
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| MISSILE 





Figure 37. Pre-Computation for Calculating BRGC 
Bis First Maneuver 
The next object is to calculate the first course to steer 
to get out of the radar resolution cell after having fired the 


chaff. This situation is shown in Figure 38. 


The choice to select one of the course options is based 
on the relative movement of the ship along the vertical axis, i.e., 
XRP direction which is oriented along the missile line of flight. 
The following algorithm in the program enables us to decide what 


course option to take: 
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If YRP(N) > O then CSE2 or the first leg is into the 






missile. 
If YRP(N) < O then CSE2 is away from the missile. 
MISSILE 
CSEOPT(1) 
CSEOPT(2) 
Figure 38. Course to Steer 
4. Second Maneuver 


Once the ship is out of the radar resolution cell ve 
calculate the next course to steer to maximize miss distance. This 
is initiated when one of the following conditions occur: 

a) YRP(N) > RG 
b) XRP(N) > AG 

In either case, the object is to ascertain whether the 
ship is out of the radar resolution cell, implying that break lock 
has occurred. However, the subroutine for calculating the 
parameters for the second maneuver only calculates CSE3 and the 
turn direction. The actual implementation of affecting this course 


alteration based on the two conditions already mentioned, is 
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achieved in the main program. Since the second leg of the new 
course, CSE3, is at right angles to the first leg, CSE2, there are 
again two choices presented as shown in Figure 39. 

It has been supposed here that CSE2 is equal to CSEOPT2, 


referring to Figure 38. 


Moreover, 
CSEOPT(3) = BRGM(1) + 90 (3.112) @ 
| 
CSEOPT(4) = BRGM(1) - 90 (3.113) $ 
WIND | 
MISSILE 
| 
| 
BRGC | 
CSE2 
CSEOPT(4) CSEOPT(3) 


Figure 39. Second Leg of the Course 
The selection of CSE3 from these two options is related 
to the bearing in which the chaff is launched and is usually in its 
opposite direction. In the situation shown in Figure 39, À 
CSE3 = CSEOPT(4) 
The direction to turn is calculated in a manner similar 


to the one described in Section 3. 
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IV. DESCRIPTION OF PROGRAM 


The simulation program is written in FORTRAN-77 (Microsoft 
version 3.2). In order to show the run of play at various times, 
the program uses a compatible Plot-88 (1984) graphics package. 
A. INITIAL CONDITIONS 

The initial conditions are set up through an interactive 
exchange with the operator in order to input all the important 
variables. The following variables are used as inputs and they have 
been assigned values that are considered reasonably standard but 
they can be altered through a Simple modification in the program. 

TABLE 4 


INPUT VARIABLES WITH TYPICAL VALUES 


Variable Description Value 
CRAD Radius of Chaff 50 yards 
SL Ship Length 200 mtrs 
SB Ship Breadth 20 meter 
Max Ship Velocity 25 knots 
Missile Nav Constant 3 
Initial Ship Vel 20 knots 


Missile Velocity 600 knot 


Turn Rate of the Ship 5/s 


Ship’s Acceleration 0.5 kts/s 
Missile Acceleration 30 yards 


Range of Chaff 200 yards 
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A scaling factor is included in the input variables to enable 
actual ranges, velocities and accelerations to be depicted on the 
screen of the monitor. The screen center (4.0,3.0) is used as a 
reference with (0.0,0.0) being the lower left hand corner and 


(8.0,6.0) the upper right hand corner of the screen. 


B. RUN TIME 

The program is incremented in steps of one second intervals 
until the missile gets to within 500 yards from the ship. 
Thereafter, to show details of action, we use 0.1 second intervals. 


The total time to go is found by the following relationship: 


TTG = Am (4.1) 
where 
RANGE = Initial Range of the Missile 


The total program step intervals including 0.1 second steps are 


given by: 
ITTG = Integer { TTG + (40) (500 | (4.2) 
E: MAIN BODY OF THE PROGRAM 


The main body of the program iterates through an array with 
a maximum of 60 iterations allowed. The values of missile, chaff 
and ship parameters are updated along with other related parameters 


in each time step. The related parameters include gates, angles, 
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equivalent target parameters and radar cross sections, etc. The 
missile radar seeker is turned off when its range from the 
equivalent target is 250 vyards resulting in zero missile 
acceleration and from thereon, the missile motion is essentially 
ballistic. When the missile has flown through the chaff cloud, its 
range starts increasing again and the main loop is completed on 
that iteration. If the missile passes within 50 yards of the ship, 
it is considered to have hit the ship. The implementation of 
alteration of courses to maneuver the ship is also achieved in the 
main body of the program using the decision algorithms mentioned 


in Section ' L’. 


D. USER INTERACTION AND OUTPUT 

After having initialized all the variables, the user is 
presented with two options; either to let the program calculate the 
control parameters by resorting to auto mode or to manually inject 
these parameters. The parameters consist of chaff firing bearing, 
both legs of the courses to be steered to maneuver the ship and the 
turn directions. Subsequent to this stage, the user is presented 
with a picture comprising of initial conditions as shown in Figure 
40. The essential parameters are listed on the picture to its 
right and at the bottom. Thereafter, the user is presented with 
another menu of the following options: 

T. Option One 

This option draws the missile, ship and chaff paths and 


gives the miss distance at the bottom as shown in Figure 41. 
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Initial Conditions 


Figure 40. 
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Figure 41. Paths Traversed by Missile, Ship & Chaff 
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TARDO 


036,1 


MINIMUM RANGE TO SHIP 


2. Option Two 


This option provides a pictorial representation quite. 


similar to the initial conditions picture and we see the situation 
in the resolution cell in 5 second increments. It depicts the run 


of play showing successive positions of the missile, chaff and the 


ship. An example is shown in Figure 42. 





= 


Figure 42. Instantaneous Position of the 
Resolution Cell 
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parameters with the help of a menu. It is incorporated to provide 
flexibility to the user to make any alterations in the input 
parameters as desired. 


either individually or in totality. The menu of options is as 


3. Option Three 


Option three facilitates changes in the input scenario 


follows: 


into smaller modules. Their respective functions and broad outlines 


are presented here to give an overview to the reader. 


TABLE 5 


CHANGEABLE PARAMETERS 


MISSILE BEARING 

MISSILE RANGE 

SHIP’S HEAD 

SCREEN SCALE 

MISSILE RADAR BEAMWIDTH 

RANGE GATE 

MAXIMUM SHIP RADAR CROSS SECTION 
MINIMUM SHIP RADAR CROSS SECTION 
RADAR CROSS SECTION OF CHAFF 
TRUE WIND DIRECTION 

TRUE WIND SPEED 

ALL PARAMETERS 

RETURN TO MAIN PROGRAM 





4. Option Four 


This option terminates the progran. 


SUBROUTINES 


A number of subroutines are made use of to break the program 
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The parameters can be changed by choice 


E Auto Subroutines 
There are two auto subroutines. The function of the first 
auto subroutine is to calculate the chaff launch bearing, BRGC, the 
first leg of the ship maneuvers, CSE2 and the direction to enable 
the shorter turn from the initial course. In order to optimize miss 
distance, a fairly accurate estimate of the chaff launch bearing 
is obtained using the pre-computation method. The decision related 
to CSE2 calculation is based on the algorithm in section “L“. The 
second subroutine works out the second leg of the maneuver and the 
direction to turn away from the chaff. 
28 Ship, Rescell And Chaff Subroutines 
These subroutines calculate the coordinates of the 
corners of the resolution cell, the hexagon representing the chaff 
cloud and the polygon representing the ship. 
3. Shipcell And Chaffcell Subroutines 
These subroutines calculate the radar cross sections of 
the ship and the chaff inside the resolution cell to enable 
calculation of the centroid in the main program. 
4. Aspect Subroutine 
The aspect angle of the ship as it exits the resolution 
cell is passed on to the Shipcell subroutine for subsequent RCS 
calculations. 
5; Missile, Wind, Path And Table 1 Subroutines 
These are PLOT-88 subroutines. The first three are used 


to draw the missile line, wind line and paths traced by the 
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missile, ship and chaff. Table 1 subroutine constructs a table of 
essential parameters when Screnl & 2 subroutines are used. 
6. Scren 1 2 & 3 Subroutines 

Scren 1 subroutine draws the missile line of flight, wind 
vector and the ship polygon based on the initial conditions input 
by the operator. Scren 2 subroutine is used to draw the same 
picture but at 5 second intervals to offer a look at the situation 
dynamically. Scren 3 subroutine enables the operator to change 
desired parameters and updates the memory with their latest 


version. 
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V.  RESULTS/RECOMMENDATIONS 


A- RESULTS 

The random variability of the parameters associated with the 
ship, missile and the chaff in a given battle scenario can result 
in a very large number of possible situations at sea. Some 
combinations of the parameters mentioned can produce a more 
favorable situation for the ship than others to thwart the anti- 
ship missile. In some cases, a single chaff round deployed by 
itself would have reduced effectiveness or might not be effective 
at all. A strong wind blowing from the direction of the missile or 
from the opposite side will move the chaff via the shortest path 
out of the resolution cell. In such cases, the chaff will exit the 
cell before the ship can do so; thus foiling the effort to break 
lock. Firing the chaff into the predicted range gate would avoid 
this problem. It is not considered appropriate to discuss results 
produced by the algorithm used in this work for all the possible 
cases. Some random situations have, therefore, been selected to 
give the reader an idea of the working of the model developed. 

1. Situation 1 


The salient input variables for this situation are as 


follows: 
e Missile Bearing - 000 Degrees 
e Missile Range - 10,000 Yards 
e Ship’s Course - 225 Degrees 
e Wind Direction - 330 Degrees 
e Wind Speed = 15 Knots 
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The initial picture is shown in Figure 43. The bearing to launch 
as calculated by the program when selected in automatic mode is 
090. This has been worked out by the pre-computation technique. If 
the manual mode is selected instead, the user is presented the same 
picture without the launch bearing and asked to input the 
following: 

Bearing to launch the chaff 

First course to steer to exit the range gate 

Direction to turn to first course 

e Second course to steer to get away from the missile line of 


fire 
e Direction to turn to second course 


In the automatic mode, however, the program works out the 
solution itself and then the result of the simulation can be viewed 
in two ways. The first graphic display shows the paths traced by 
the ship, the missile and the chaff. This is shown in Figure 44. 
In the case shown, the missile has missed the ship by 882.9 yards. 
The ship first maneuvered to steer southwards to break lock and 
then westwards to maximize miss distance. The second mode of 
presentation displays the situation every 5 seconds until the 
missile hits one of the two possible targets. The situation after 


15 seconds is shown in Figure 45. 
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Initial Conditions(Situation 1) 


Figure 43. 
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Figure 44. Paths display (Situation 1) 
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Figure 45. Instantaneous Positions(Situation 1) 
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2. Situation 

The input variables 

ə Missile Bearing 

e Missile Range 

e Ship’s Course 

e Wind Direction 

e Wind Speed 
Figure 46 shows the 


chaff is 010. The 


2 


are as follows: 


280 Degrees 
10,000 Yards 
150 Degrees 
220 Degrees 
15 Knots 


initial conditions. The bearing to launch the 


paths traversed and the situation after 15 


seconds are shown in Figures 47 and 48 respectively. 
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Initial Conditions(Situation 2) 


Figure 46. 
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Figure 47. Paths display (Situation 2) 
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Figure 48. Instantaneous Positions(Situation 2) 
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B. RECOMMENDATIONS 

The deployment of chaff in seduction mode poses a complicated 
problem in terms of the randomness of those variables that are not 
controllable by the ship. They include the missile“s bearing, the 
beam width of its radar, the ship“s course at the time of initial 
detection and the wind vector. Combined with the initial detection 
range of the missile these factors determine the size of the radar 
resolution cell and as a result, the time spent by the chaff round 
inside the cell. This is also the time duration available to the 
ship to vacate the cell and is, therefore, critical in solving the 
problem. 

The variables that can be controlled by the ship are the chaff 
range and bearing, ship's course and speed, its turn rates and 
acceleration. The last three have a limited contribution in this 
context. The control of the ship's heading and the bearing to fire 
the chaff is of the utmost significance. In order to increase the 
effectiveness of the control variables and to make the algorithm 
developed more adaptable, some recommendations are made in the 
ensuing paragraphs. 

ão: Multiple Chaff Firing 

In order to tackle the cases where the chaff exits the 
radar resolution cell before the ship, multiple chaff rounds must 
be launched. The second chaff should be reseeded inside the cell 
prior to the exit of the first. Its initial placement should depend 
on how far the ship is from the cell bounds when the first round 


is about to exit or how much time it requires at the point to leave 
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the cell keeping the bloom time delay in mind. Thereafter, 
subsequent chaff rounds may be fired in the same way as required 
until the ship breaks lock. 

2: Chaff Walk off 

This is another philosophy involving multiple chaff 
launches. If multiple chaff rounds are fired on the same bearing 
but at increasing ranges, rapid break lock could be achieved and 
then the ship would need to carry out only one maneuver to increase 
miss distance. The chaff rounds can be fired simultaneously or in 
quick succession so that they are waiting in anticipation and the 
centroid shifts in a discrete manner. To explain this, let us study 
an example shown in Figure 49. 

With the missile approaching from the north, the first 
chaff is placed at the cell's angle gate boundary on a bearing 
orthogonal to missile direction. The bearing can be offset to cater 
for the bloom time delay. If the chaff radar cross section is equal 
to that of the ship, the centroid will be established at 
approximately the mid-point of the two, called X1. The resolution 
cell now moves to the right and is called Cell ’B’. At this point 
it includes the ship, Chaff 1 and Chaff 2. The centroid now shifts 
to X2 and the cell to Cell ’C’. The third chaff round displaces 
the cell position to Cell ’D’ and the centroid to X3. It can be 
seen that the ship is now very close to the left hand angle gate 
boundary and can steer out of the cell to achieve a rapid break and 
then increase the miss distance. We feel that this technique can 


enhance the effectiveness of chaff in anti-ship missile defense. 
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Figure 49. Chaff Walk Off 


a). Multiple Solutions 
In order to achieve optimum results, the program can be 
used in manual mode to inject different control inputs. The program 
run time is quite short and a slight modification to compare 
results of both automatic and manual modes can recommend the best 
solution to be adopted in a particular situation. 
4. Battle Group Operation 
When operating in a battle group,  non-judicious 


deployment of chaff can result in seducing the chaff onto a 
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friendly unit. It is, therefore, imperative to introduce safe. 


firing bearings in the algorithm to ensure safety of ships in 


company. 
5: Fixed Chaff Launcher 
The program needs to be adapted to include units equipped 
with fixed chaff launchers to optimize firing bearings while 
minimizing maneuvers by the ship. The complexity of the problem is, 
of course, increased and the ship’s actions are further limited. 
If there are limitations of CIWS firing bearings, they can also be 
included to allow for an integrated response. 
6. Integration with Active ECM 
If the ship is equipped with an active ECM capability, 
integration of chaff deployment with the active ECM response to an 
anti-ship missile threat can be a decisive factor. The break lock 
may first be achieved with active ECM and chaff can then be placed 
between the ship and the false target to deceive the missile and 
also to increase miss distance. Active decoys such as Nulka and 
Winnin employ similar techniques. 
7. Stochastic Modelling 
Due to the randomness of the variables involved in the 
anti-ship defense problem, stochastic modelling may be employed to 
achieve optimum results. In this paper, wind direction and speed 
are assumed to be constant. This is usually not true in real sea 
conditions. Noise in the wind vector can also be introduced through 
an appropriate probability distribution to depict more realistic 


environmental conditions. 
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8. Radar Cross Section Calculations 
The mathematical modelling to calculate radar cross 
sections of both the ship and the chaff, use some approximations 
when a simultaneous exit from the angle and range gates takes 
place. Also, in case of simultaneous spill-over from opposing sides 
of the cell through the range gate, the partial radar cross section 
is calculated assuming that the spill-over is from one side only. 
The cross section which lies outside the cell on the opposite side 
is thus neglected. The reader should be aware of these 
approximations and make suitable changes if so desired. 
9. Programming Modifications 
The simulation program is written in FORTRAN which has 
some inherent limitations. It is recommended that a more powerful 
and flexible language be used to improve upon the program 
efficiency and to incorporate the changes suggested more 
effectively. The dimension figure of the arrays can be increased 
to include larger initial detection ranges which may be applicable. 
Parameters such as chaff radius, ship's length, turning rates and 
acceleration/deceleration tables may be included as input variables 


to adapt the program to given chaff and ship charateristics. 


C: CONCLUSIONS 

Resurgence of decoys in recent years has highlighted their 
importance in the area of anti-ship missile defense. Although the 
missile radars are becoming smarter by including advanced ECCM 


features, chaff is considered to be a viable defensive technique 
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even if deployed in a stand alone manner such as centroid mode. 
This is especially valid for the smaller navies who have limited! 


resources. Improvement in chaff deployment algorithms and, when 






possible, integrated operations with active ECM and CIWS assets can. 


make chaff a potent defensive measure. It is simple to use and. 
relatively inexpensive among the known type of decoys. Its. 
intelligent use can make it effective and dependable. Although wars — 
have not been won by a single factor alone, chaff can certainly 
become one of the important factors in a limited conventional 


conflict. 
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APPENDIX. SIMULATION PROGRAM 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


$ DEPLOYMENT OF CHAFF IN CENTROID * 
* MODE AGAINST ANTI SHIP MISSILES * 
* USING A * 
* VARIABLE AZIMUTH AND ELEVATION * 
* LAUNCHER * 
* 24 AUG 89 SIDDIQ.FOR * 
. 
* 


M. SIDDIO 
e de de de de e de de de de de ER de de de de ke ke ke ke ke k k ke ke k k k k k kk k kkk kkk 


THE PROGRAM HAS NO PRINTER/ PLOTTER OPTIONS 

ITS OUTPUT GOES TO THE MONITOR ONLY 

IT HAS AN AUTOMATIC OR MANUAL CHAFF LAUNCHING OPTION 
IT CREATES NO DATA FILE 


VARIABLES DECLARATION 


COMMON/TIME/ITTG,TTG 

COMMON/TME/T 

COMMON/TRACK/ XM(60),YM(60),XS(60),YS(60),XC(60),YC(60) 
COMMON/WD/GW, VW 

COMMON/MISS/RG,THETA,AG,PI 
COMMON/ABCD/BRGM(60) ,R(60) ,SC, VS(60) , CSE(60) , BRGMT (60), 


*BRGTM(60) , RMT (60) 


COMMON/EQ/ XET(60) , YET(60) 
COMMON/SHCL/AGS , RGS, BRGMS (60) , SAA, PARTS, SAAN (60) 
COMMON/CFCL/PARTC, CRAD, RGC, AGC 

COMMON/SHP/SL, SB 
COMMON/RC/RCSS,RCSSP (60) ,RCSC,RCSCP (60) 
COMMON/MINN/RMIN 

COMMON/CHF/BRGC, ISHOOT 

COMMON/TURN/TR, CSEOFF 

COMMON/ FG/IFLAGI1 (60), IFLAG2 (60) 
COMMON/PRN/IOPORT, MODEL 


COMMON/DOTS/XMD(60),YMD(60),XMDD(60),YMDD(60),XSD(60),YSD(60), 


*XSDD(60),YSDD(60) 


COMMON/REL/XR (60) , YR(60) , XRD(60) , YRD(60) , XRP(60) , YRP(60) , XRPD(60), 


*YRPD(60) 


REAL AMM,NAV,T 
REAL CDIST 

REAL TTG,GT,X,Y 
REAL SIG,SIGD,AM,GM 
IOPORT = 99 

MODEL = 99 
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C LIST OF VARIABLES 

C 

C MISSILE: 

E p ow om > > 

C BRGM BEARING OF MISSILE FROM SHIP(WRT NORTH) 
C BRGMS BEARING OF SHIP FROM MISSILE(WRT REFERENCE) 
C BRGMT BEARING OF EOUIVALENT TARGET FROM 
MISSILE (REF) 

C BRGMC BEARING OF CHAFF FROM MISSILE(REF) 

C R/R(N) INSTANTANEOUS MISSILE RANGE FROM SHIP 
C RANGE INITIAL MISSILE RANGE 

C R(1)/RMIN SCALED INITIAL MISSILE RANGE 

e RMT (N) MISSILE RANGE FROM EQUIVALENT TARGET 
E RMC MISSILE RANGE FROM CHAFF 

E VMM MISSILE VELOCITY 

C VM SCALED MISSILE VELOCITY 

C GM MISSILE DIRECTION 

C AM MISSILE ACCELERATION 

C AMM SCALED MISSILE ACCELERATION (30/SC) 

© XM X-COORD OF MISSILE POSITION 

C YM Y-COORD OF MISSILE POSITION 

€ XMD X-COORD OF MISSILE VELOCITV 

C YMD Y-COORD OF MISSILE VELOCITY 

E XMDD X-COORD OF MISSILE ACCELERATION 

C YMD Y-COORD OF MISSILE ACCELERATION 

C SHIP 

C => «=> => «=> «> 

C SL SHIP LENGTH 

C SB SHIP BREADTH 

€ CSE COURSE 

C CSE2 FIRST LEG OF NEW COURSE TO STEER 

€ CSE3 SECOND LEG OF NEW COURSE TO STEER 

C CSEOPT COURSE OPTION FOR FIRST/ SECOND LEG 

C CSEDIF1 DIFFERENCE BETWEEN CSE1 AND CSE2 

C TOPT TIME TO TURN TO CSEOPT 

C as ANGLE BETWEEN REF AND SHIP’S HEAD 

C VS SCALED SHIP VELOCITY 

E AS (VR) SHIP ACCELERATION 

E TR TURN RATE OF SHIP 

C U(ALSO U1,U2,U3)TURN DIRECTION OF SHIP 

C SAA ASPECT ANGLE OF SHIP(BETWEEN CSE & BETA) 
C XS X-COORD OF SHIP POSITION 

E YS Y-COORD OF SHIP POSITION 

C XSD X-COORD OF SHIP VELOCITY 

C YSD Y-COORD OF SHIP VELOCITY 

C XSDD X-COORD OF SHIP ACCELERATION 

© YSSD Y-COORD OF SHIP ACCELERATION 

€ RCSS RADAR CROSS SECTION OF SHIP 

C RCSSMAX MAXIMUM RADAR CROSS SECTION OF SHIP 

C RCSSMIN MINIMUM RADAR CROSS SECTION OF SHIP 

C RCSSP(1) INITIAL RADAR CROSS SECTION OF SHIP 

C RCSSP(N) INSTANTANEOUS RADAR CROSS SECTION OF SHIP 
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© BETA ANGLE BETWEEN REF USED FOR MICROSOFT 
FORTRAN 6 
e BRGM 
C RESOLUTION CELL: 
ČC -a 
€ RGATE RANGE GATE 
C RG SCALED RANGE GATE 
© AG ANGLE GATE 
€ RGS DISTANCE BETWEEN RESCELL CENTER € SHIP 
CENTER 
C ALONG THE RANGE GATE 
€ RGC 
RR === e E EE COAEF======= 
€ AGS DISTANCE BETWEEN RESCELL CENTER € SHIP 
CENTER 
C ALONG ANGLE GATE 
C AGC 
A == === === > == === == == GRE == 
G S2r/S2a DISTANCE BETWEEN RG/AG BOUNDARY & SHIP 
CENTER 
C Rir/RI1a 
eene e EE EE e CHAFF------- 
Sir/Sla PORTION OF HALF OF THE LONGITUDINAL AXIS 
INSIDE THE CELL 
PARTS PORTION OF THE SHIP INSIDE THE RESCELL 
MT TEMMA l... EA I oe 
CHAFF: 
BRGC BEARING TO LAUNCH CHAFF 
CDIST CHAFF LAUNCH DISTANCE 
CRAD CHAFF RADIUS 
VC CHAFF SPEED 
GC DIRECTION OF MOTION OF CHAFF 
BRGC1&2 BEARINGS PERPENDICULAR TO BRGM TO FIRE 
HAFF 
RCSC RADAR CROSS SECTION OF CHAFF 
RCSSC (1) INITIAL RADAR CROSS SECTION OF CHAFF (=0) 
RCSSC (N) INSTANTANEOUS RADAR CROSS SECTION OF CHAFF 
XC X-COORD OF CHAFF POSITION 
YC Y-COORD OF CHAFF POSITION 


EQUIVALENT TARGET: 


ARGET 


009000000000000000000000000000 


WIND: 


BEARING OF MISSILE FROM EQUIVALENT TARGET 
DIRECTION OF MOTION OF EQUIVALENT TARGET 


EQUIVALENT TARGET ANGLE 

X-COORD OF EQUIVALENT TARGET POSITION 
Y-COORD OF EQUIVALENT TARGET POSITION 
RATE OF CHANGE OR VELOCITY OF EQUIVALENT 


POSITION 
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YRP 


XRPD(1) 
YRPD(1) 
YRPD(N) 
ELOCITY 


Q<4000000000000000000 


XRPD(N) 
VELOCITY 

c XRP1 
DIRECTION 

© YKP1 
DIRECTION 

C XRP2 
DIRECTION 

e YRE2 
DIRECTION 
XRPD1 
YRPD1 
XRPD2 
YRPD2 
ARPDD 
YRPDD 


TTG 
ITTG 
XCENT 
YCENT 
NAV 
SIG 
IGHT (CLOCKWISE) 
SIGD 


PAANAQAAQNQAANAANAIANAANAIANNA 


O 


NAV = 3.0 
T=0 

H = 1. 
XCENT 


I o 
b 
O 


MISCELLANEOUS: 


WIND DIRECTION(FROM WHERE IT COMES FROM) 


TRUE WIND SPEED 
SCALED WIND SPEED 


TRUE WIND DIRECTION(TO WHERE IT GOES) 


RELATIVE COORDINATES 


X-COORD OF RELATIVE POSITION OF SHIP 
X-COORD OF RELATIVE POSITION OF SHIP 
X-COORD OF RELATIVE POSITION OF SHIP 
WITH THE RELATIVE COORDINATE 

SYSTEM ROTATED TO MISSILE BEARING 
X-COORD OF RELATIVE POSITION OF SHIP 
WITH THE RELATIVE COORDINATE 

SYSTEM ROTATED TO MISSILE BEARING 


INITIAL X-COORD OF RELATIVE SHIP'S VELOCITY 
INITIAL Y-COORD OF RELATIVE SHIP'S VELOCITY 


INSTANTANEOUS Y-COORD OF REL SHIP’S 


INSTANTANEOUS X-COORD OF REL SHIP’S 


DISTANCE TRAVELLED BY SHIP IN POS XRP 


DISTANCE TRAVELLED BY SHIP IN POS YRP 


DISTANCE TRAVELLED BY SHIP IN NEG XRP 


DISTANCE TRAVELLED BY SHIP IN NEG YRP 


X-COORD OF RELATIVE 
Y-COORD OF RELATIVE 


SCREEN SCALE 
TIME TO GO 
PROGRAM RUN TIME 


RELATIVE VELOCITY OF SHIP IN XRP1 DIRECTION 
RELATIVE VELOCITY OF SHIP IN YRP1 DIRECTION 
RELATIVE VELOCITY OF SHIP IN XRP2 DIRECTION 
RELATIVE VELOCITY OF SHIP IN XRP2 DIRECTION 


ACCELERATION OF SHIP 
ACCELERATION OF SHIP 


X-COORD OF SCREEN CENTER 
Y-COORD OF SCREEN CENTER 
PROPORTIONAL NAVIGATION CONSTANT 
ANGLE BETWEEN REF AND LINE OF 


RATE OF CHANGE OF LOS 


PI = 4.0 * ATAN(1.0) 
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YCENT = 3.0 


SL = 200.0 
SB = 20.0 
IFLAG1(1) = 0 
IFLAG2(1) = 0 
IFLAG3 = 0 
C 
WRITE(*,*) /INPUT MISSILE BEARING (000-359 DEGREES TRUE) “ 
READ(*,*) BRGM(1) 
WRITE(*,*) “INPUT INITIAL MISSILE RANGE (< 15000 YARDS) ’ 
READ(*,*) RANGE 
WRITE(*,*) ‘INPUT SHIP HEADING (000-359 DEGREES TRUE) ’ 
READ(*,*) CSE(1) 
C MISSILE PARAMETERS 
C wem em em am an an an an an an am an am am am am amw amw 
WRITE(*,*) “INPUT MISSILE RADAR BEAMWIDTH (+/- 1.0 TO +/-5.0 
*DEGREES) ' 
READ(*,*) THETA 
WRITE(*,*) ‘INPUT THE RANGE GATE (+/-10 TO +/-100 YARDS) ’ 
READ(*,*) RGATE 
C RADAR CROSS SECTION 
© — mm em em e em em em em em em em em e ep em ep em 
G SHIP RCS ASSUMED MAX RIGHT AHEAD/ASTERN & ABEAM AND MIN AT 
BOWS/QTRS 


WRITE(*,*)’INPUT RADAR CROSS SECTION OF SHIP & CHAFF’ 
WRITE(*,*)’ FOR CHAFF TO BE EFFECTIVE THERE MUST BE AT LEAST’ 
WRITE(*,*)’ A 1:2 RCS RATIO (CHAFF:SHIP) ’ 
WRITE (*,*)" 7 
WRITE(*,*) “INPUT MAX RADAR CROSS SECTION OF SHIP ’ 
WRITE(*,*)/10000 - 100000 SQUARE METERS’ 
WRITE(*,*)” 7 
READ(*,*) RCSSMAX 
WRITE(*,*) “INPUT MINIMUM RADAR CROSS SECTION OF SHIP’ 
WRITE(*,*)/5000 - 50000 SQUARE METERS’ 
WRITE(*,*)”  / 
READ(*,*) RCSSMIN 
WRITE (*,*) "INPUT RADAR CROSS SECTION OF CHAFF’ 
WRITE(*,*)'5000 - 20000 SQUARE METERS ’ 
WRITE(*,*)” 7 
READ(*,*) RCSC 

C WIND PARAMETERS 


C em zm mm wm mm mm mm zm mm mm zm mm zm zm em 
WRITE(*,*) /INPUT WIND DIRECTION (000-359 DEGREES TRUE) “ 
READ(*,*) GW 
WRITE(*,*) “INPUT TRUE WIND SPEED (0.0-50.0 KNOTS) ” 
READ(*,*) VW 
WRITE(*,*) “INPUT SCREEN SCALE (YDS:1 INCH) ” 
WRITE(*,*)’ RECOMMEND 200 YARDS TO THE INCH’ 
READ(*,*) SC 

20 CONTINUE 

ee SCALING 

e —— a 
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VS = SHIP'S SCALED VELOCITY CONVERTED TO YDS/SEC 
VS(1) = 20./(SC*1.8) 


RG = SCALED RANGE GATE 
RG = RGATE/SC 

VM = SCALED MISSILE VELOCITY 
VMM = 600. 


VM = VMM/(SC*1.8) 

VW1 = SCALED WIND VELOCITY 
VW1 = VW/(SC*1.8) 

SCALE THE RANGE FOR PLOTTING 


R(1) = RANGE/SC 

RMIN = R(1) 

TR IS THE TURN RATE IN DEGREES/SEC 
TR = 5.0 


VR IS THE ACCELERATION IN KNOTS/SEC*SC 
VR = 0.5/(SC*1.8) 


RCSS=RCSSMIN+ABS ( (RCSSMAX-RCSSMIN) *COS (2* (CSE(1) -BRGM(1) ) *PI/180) ) 


aa 


aa 


0090 


RCSSP(1) = RCSS 
RCSCP(1) = 0.0 


TTG = TIME TO GO 


HI 
wel 
Q 

Il 


RANGE*1.8/VMM 
ITTG = IFIX(TTG +10*500./(VM*SC)) 


MISSILE ACCELERATION 


AMM=30./SC 


SIG REF. TO HORIZONTAL (EAST) 


SIG = (270. - BRGM(1))*PI/180. 
IF (BRGM(1).GT.180.) THEN 
GM = BRGM(1) - 180. 
ELSE 
GM = BRGM(1) + 180. 
END IF 


XXXXXXINITIAL CONDITIONS *žíží* 
MISSILE I. C. 


XM(1) = XCENT + R(1)*SIN(BRGM(1)*PI/180.) 
YM(1) = YCENT + R(1)*COS(BRGM(1)*PI/180.) 


XMD(1) = VM*SIN(GM*PI/180.) 
YMD(1) = VM*COS(GM*PI/180. ) 
SHIP I C. 


XS(1) = XCENT 
YS(1) = YCENT 
XSD(1) = VS(1) * SIN(CSE(1)*PI/180.) 


YSD(1) = VS(1) * COS(CSE(1)*PI/180.) 
XSDD(1) = 0.0 
YSDD(1) = 0.0 


GS = ATAN2(YSD(1),XSD(1)) 


FIND THE NEW COURSE TO STEER AND BEARING TO LAUNCH CHAFF 
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—< E 


30 WRITE (*,*)” INPUT (1) FOR AUTOMATIC CHAFF LAUNCH ' 
WRITE (*,*)” AND NEW SHIP COURSE ’ 
WRITE(*,*)” 7 
WRITE (*,*)” INPUT (2) FOR MANUAL CHAFF LAUNCH ’ 
WRITE (*,*)' AND NEW SHIP COURSE ' 
READ(*,*) ISHOOT 
IF (ISHOOT.EQ.1) THEN 

CALL AUTO1 (CSE2,U2) 


U = U2 
CALL AUTO2 (CSE3,U3,CSE2) 
Ul = U3 
© DRAW THE INITAL MISSILE RESOLUTION CELL 


WRITE(*,*)/**** MISSILE DETECTION **** / 
RMT(1) = R(1) 
BRGTM(1) = BRGM(1) 
XET(1) = XS(1) 
YET(1) = YS(1) 
CALL SCREN1 
ELSE IF (ISHOOT.EQ.2) THEN 
WRITE(*,*)'**** MISSILE DETECTION *x*xx f 
RMT(1) = R(1) 
BRGTM(1) = BRGM(1) 
XET(1) = XS(1) 
YET(1) = YS(1) 
CALL SCREN1 
WRITE (*,*) “INPUT TRUE BEARING TO LAUNCH CHAFF’ 
READ(*,*) BRGC 
WRITE(*,*)' INPUT FIRST LEG OF NEW SHIP COURSE TO STEER ’ 
READ(*,*) CSE2 
40 WRITE(*,*)” INPUT (1) TO TURN THE SHIP RIGHT ’ 
WRITE (*,*)” INPUT (-1) TO TURN THE SHIP LEFT ’ 
READ (*,*) U2 
U = U2 
IF ((U2.NE.1.0).AND.(U2.NE.-1.0)) THEN 
GO TO 40 
END IF 
WRITE (*,*) "INPUT SECOND LEG OF NEW SHIP COURSE TO STEER ’ 
READ (*,*) CSE3 
50 WRITE (*,*)” INPUT (1) TO TURN THE SHIP RIGHT ’ 
WRITE (*,*)” INPUT (-1) TO TURN THE SHIP LEFT “ 
READ (*,*) U3 
Ul = U3 
IF ((U3.NE.1.0).AND. (U3.NE.-1.0)) THEN 
GO TO 50 
END IF 
ELSE 
GO TO 30 
END IF 
C CHAFF INITIAL CONDITIONS 
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C © CDIST = DISTANCE FM SHIP CHAFF IS LAUNCHED 
C CRAD = CHAFF RADIUS 
C CHAFF IS FIRED TO THE EDGE OF THE CELL BUT THE DISTANCE IS 
REDUCED 
C BY ONE CHAFF RADIUS TO CATER FOR THE SHIP'S MOVEMENT DURING THE 
C BLOOM TIME DELAY 
C VC = CHAFF SPEED IN YDS/SEC 
CRAD = 100./SC 
CDIST = (R(1)*TAN(THETA * PI/180)) - CRAD 
VC = VW / (SC*1.8) 
GC = GW + 180. 
IF (GC.GT.360) THEN 
GC = GC - 360 
ENDIF 
XC(1) 
YC(1) 


XS (1) 
YS(1) 


ann 


kkkkkkkkkkkk MAIN  LOOP *kkkkkxxkkkkxk xxx x 


DO 100 N = 1, ITTG 
IF (IFLAG1(N-1).EQ.1) GOTO 110 


ann 


IFLAG1(N) = O 
IFLAG2(N) = O 
C RANGES AND ANGLES (WRT HORIZONTAL REFERENCE 


R(N) = SORT((XS(N)-XM(N))**2+(YS(N)-YM(N))**2) 

RCM = SORT((XC(N)-XM(N))**2+(YC(N)-YM(N))**2) 

RMT(N) = ((XM(N)-XET(N))**2+(YM(N)-YET(N))**2)**0.5 
IF(((XET(N)-XM(N)).EQ.0.).AND. ((YET(N)-YM(N)).EQ.0.)) THEN 


BRGMT(N) = 0.0 
ELSE 

BRGMT(N) = ATAN2(YET(N)-YM(N),XET(N)-XM(N)) 
END IF 


IF(((XET(N)-XM(N)).EQ.0.).AND.((YET(N)-YM(N)).EQ.0.))THEN 
BRGTM(N) = 0.0 

ELSE 
BRGTM(N) = (180./PI) *ATAN2 (XM(N) -XET(N) , YM(N) -YET(N) ) 

END IF 

IF (((XS(N) -XM(N)).EQ.0.).AND. ((YS(N)-YM(N)).EQ.0.)) THEN 
BRGMS(N) = 0.0 

ELSE 
BRGMS (N) = ATAN2(YS(N)-YM(N),XS(N)-XM(N)) 

END IF 

IF(((XC(N)-XM(N)).EQ.0.).AND. ((YC(N)-YM(N)).EQ.0.)) THEN 
BRGMC = 0.0 

ELSE 
BRGMC = ATAN2 (YC(N)-YM(N),XC(N)-XM(N)) 

END IF 

C RANGE AND ANGLE GATES 
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AGS = R(N) * SIN(ABS(BRGMS (N) -BRGMT (N)) ) 

AGC = RCM * SIN(ABS(BRGMC-BRGMT (N))) 

RGS = ABS(R(N) * COS(ABS(BRGMS (N) -BRGMT (N) ) ) -RMT(N)) 
RGC = ABS(RCM * COS (ABS (BRGMC-BRGMT (N) ) ) -RMT (N) ) 


AG = RMT(N) * TAN(THETA*PI/180.) 


RCSS=RCSSMIN+ABS ( (RCSSMAX-RCSSMIN) *COS (2* (CSE(N) -BRGM(N)) 
* *PI/180)) 
RCSSP(N) = RCSS 
IF ((RGS.GT.(RG-SL/(SC*2.0))).OR. 
* (AGS.GT.(AG-SL/(SC*2.0)))) THEN 
CALL SHPCLL(N) 
RCSSP(N) = RCSS * PARTS 
END IF 
C CHAFF RCS IS ASSUMED TO BE ZERO UNTIL FULL BLOOM 


IF (N.GT.5) THEN 
RCSCP(N) = RCSC 
IF ((RGC.GT. (RG-CRAD) ).OR. (AGC.GT. (AG-CRAD) )) THEN 
CALL CHFCLL 
RCSCP(N) = RCSC * PARTC 
END IF 
ELSE 
RCSCP(N) = 0.0 
ENDIF 
C EQUIVALENT TARGET 


IF ((RCSSP(N).EQ.0.0) .AND. (RCSCP(N) .EQ.0.0)) THEN 
XET(N+1) = XET(N) 
YET(N+1) = YET(N) 
ELSE 
XET (N+1) 
* RCSCP(N)) 
YET(N+1) = (RCSSP(N)*YS(N)+RCSCP(N)*YC(N))/(RCSSP(N)+ 
* RCSCP(N)) 
END IF 
C ANGLE BETWEEN HORIZONTAL REFERENCE AND DIRECTION OF MOVEMENT 
OF 
C EQUIVALENT TARGET 
| IF (N.GE.2) THEN 


(RCSSP(N) *XS(N)+RCSCP(N) *XC(N))/(RCSSP(N)+ 


IF(((XET(N)-XET(N-1)).EQ.0.).AND. ((YET(N)-YET(N-1)).EQ.0.)) 
k 


THEN 
GET = 0.0 
ELSE 
GET = ATAN2(YET(N)-YET(N-1),XET(N)-XET(N-1)) 
END IF 
C EQUIVALENT TARGET VELOCITY 


VET = SORT((XET(N)-XET(N-1))**2+(YET(N)-YET(N-1))**2)/H 
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ELSE 


GET = GS 
VET = VS(N) 
END IF 


C MISSILE NAVIGATION DYNAMICS 


AR 


IF(((XET(N)-XM(N)).EQ.O.).AND. ((XET(N)-XM(N)).EQ.O.))THEN 
SIG = 0.0 
ELSE 
SIG = ATAN2((YET(N)-YM(N)), (XET(N)-XM(N))) 
END IF 
IF (RMT(N).EQ.0.0) THEN 
SIGD = 0.0 
ELSE 
SIGD = (VET*SIN(GET-SIG) -VM*SIN(GM-SIG))/RMT(N) 
END IF 
AM=SIGD*VM*NAV 
IF(AM.GT.AMM) THEN 
AM=AMM 
ELSE IF(AM.LT.-AMM) THEN 
AM = -AMM 
END IF 


TO TURN OFF MISSILE SEEKER 


IF (RMT(N).LT.(250./SC)) THEN 


AM = 0.0 
IF (RMT(N-1).LT.RMT(N)) THEN 
ITTG = N 
IFLAG1(N) = 1 
END IF 
END IF 


E MISSILE EQUATIONS 


XMDD(N)= -AM*SIN(SIG) 
YMDD(N)= AM*COS (SIG) 
IF((YMD(N).EQ.0.).AND. (XMD(N).EQ.0.)) THEN 


GM = 0.0 
ELSE 

GM = ATAN2(YMD(N) , XMD(N) ) 
END IF 


XM(N+1)= XM(N) + XMD(N)*H + 0.5*XMDD(N) *H**2 
YM(N+1)= YM(N) + YMD(N)*H + 0.5*YMDD(N) *H**2 
XMD(N+1)= VM*COS(GM) + XMDD(N)*H 
YMD(N+1)=  VM*SIN(GM) + YMDD(N)*H 


© SHIP EQUATIONS 


IF (BRGM(N).LT.0.0) THEN 
BRGM(N) = BRGM(N) + 360. 

END IF 

XSDD(N) = -VS(N)*(PI*TR/180.)*U1*SIN(CSE(N)*PI/180.) + 


(PI*CSE(N)/180.) 
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YSDD(N) = VS(N)*(PI*TR/180.)*U1*COS(CSE(N)*PI/180.) + 


VR*COS 


* (PI*CSE(N)/180.) 
XSD(N) = VS(N)*SIN(PI*CSE(N)/180.) 
YSD(N) = VS(N)*COS(PI*CSE(N)/180.) 
IF (VS(N).LT.(25./(1.8*SC))) THEN 
VS(N+1) = VS(N) + VR*H 


ELSE 
VS (N+1)=25./(1.8*SC) 
END IF 
XS(N+1) = XS(N) + XSD(N)*H + 0.5*XSDD(N) *H**2 
YS(N+1) = YS(N) + YSD(N)*H + 0.5*YSDD(N) *H**2 


IF(((XM(N+1)-XS(N+1)) .EQ.0.) .AND. ((YM(N+1)-YS(N+1)).EQ.0.)) 
THEN 


BRGM(N+1) = 0.0 
ELSE 
B RG M (N + 1.) = 


(180./PI) *ATAN2 (XM(N+1)-XS(N+1),YM(N+1)-YS(N+1)) 
END IF 


na na 


050000 


CHAFF DYANAMICS AND RELATIVE COORDINATES EQUATIONS 


IF (N.LE.5) THEN 
CHAFF'S ARBITRARY POSITION BEFORE IT BLOOMS 


XC(N+1) = XS(N+1) + CDIST*SIN(BRGC*PI/180.) 


YC(N+1) = YS(N+1) + CDIST*COS (BRGC*PI/180.) 
ELSE 
XC(N+1) = XC(N) + VC*SIN(GC*PI/180.) 


YC(N+1) = YC(N) + VC*COS(GC*PI/180.) 
ENDIF 
SHIFT OF THE ORIGIN OF COORDINATE SYSTEM TO THE CHAFF 
= XS(N)-XC(N) 
YR(N) = YS(N)-YC(N) 
ROTATION OF THE COORDINATE SYSTEM TO THE MISSILE DIRECTION 
XRP(N) = (XR(N)*COS(BRGM(N)*PI/180)) - (YR(N)* 
* SIN(BRGM(N)*PI/180)) 


YRP(N) = (XR(N)*SIN(BRGM(N)*PI/180)) + (YR(N)* 
* © COS(BRGM(N)*PI/180)) 
T=T+H 


ALTERATION OF COURSES 
WATCH FOR CSE(N) - CSE2/CSE3 NEVER BEING LESS THAN TR 
WHILE THE SHIP IS WITHIN THE RG BOUNDS UPDATE COURSE TO ALTER 


THE FIRST LEG UNTIL THE TURN IS COMPLETED 
IF (ABS(YRP(N)).LE.RG) THEN 
IF ((ABS(CSE(N)-CSE2)).LT.TR) THEN 
U = 0.0 
END IF 
CSE(N+1) = CSE(N) + TR*H*U 
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OG 


100 
110 


C 


WHILE THE SHIP IS WITHIN THE RG BOUNDS BUT IT HAS EXITED THE 


BOUNDS UPDATE COURSE TO ALTER TO THE FIRST LEG UNTIL THE TURN 
IS COMPLETED 
IF (ABS(XRP(N)).GT.AG) THEN 
IF ((ABS(CSE(N)-CSE3)).LT.TR) THEN 
U1 = 0.0 
ENDIF 
CSE(N+1) = CSE(N) + TR*H*U1 
ENDIF 
ELSE 
IF THE SHIP 1S OUT OF THE RG BOUNDS UPDATE COURSE TO ALTER TO 


SECOND LEG UNTIL THE TURN IS COMPLETED 
IF ((ABS(CSE(N)-CSE3)).LT.TR) THEN 
Ul = 0.0 
ENDIF 
CSE(N+1) = CSE(N) + TR*H*U1 

ENDIF 

IF (CSE(N+1).GE.360.) THEN 
CSE(N+1) = CSE(N+1) - 360.0 

END IF 

IF (CSE(N+1).LT.0.) THEN 
CSE(N+1) = 360.+CSE(N+1) 

END IF 

DETERMINE IF MISSILE HITS THE SHIP 

IF((R(N)*SC) .LT.50.0) THEN 
WRITE (*,*) ”RANGE UNDER 50 YARDS MISSILE HITS SHIP” 

END IF 

RMIN = AMIN1(RMIN,R(N)) 

IF (RMT(N).LT.(500./SC)) THEN 
IFLAG2(N) = 1 
H = 0.1 

END IF 


CONTINUE 
CONTINUE 
TO PREVENT ANY UNDEFINED VALUES IF ITTG < 60 


DO 120 I=(ITTG+1),60 


XC(I) = XC(I-1) 
XC(I) = XC(I-1) 
YC(I) = YC(I-1) 
XS(I) = XS(I-1) 
YS(I) = YS(I-1) 
XM(I) = XM(I-1) 
YM(I) = YM(I-1) 
XET(I) = XET(I-1) 
YET(I) = YET(I-1) 


94 


120 CONTINUE 


C 
C 
G 
RMIN = SC*RMIN 
Ckkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 
C PRINT & PLOT 
C em, em em em em em em em. em emm em «em 
130 WRITE(*,*)’INPUT (1) TO DRAW MISSILE & SHIP & CHAFF PATH’ 
WRITE (*,*) ’ 
WRITE(*,*) “INPUT (2) TO DRAW RESOLUTION CELL IN ’ 
WRITE (*, *) ’ FIVE SECOND INCREMENTS , 
WRITE(*,* , 
WRITE (*,*) INPUT (3) TO CHANGE SCENARIO PARAMETERS‘ 
WRITE (*, * 
WRITE (*,*) INPUT (0) TO QUIT , 


WRITE(*,*) ’ 
READ(*,*) IDRW 
IF(IDRW.EO.1) THEN 
CALL PATH 
ELSE IF(IDRW.EO.2) THEN 
CALL SCREN2 
ELSE IF(IDRW.EO.3) THEN 
CALL SCREEN3 (IFLAG3) 
IF (IFLAG3.EQ.1) THEN 
GO TO 10 
ENDIF 
IF (IFLAG3.EQ.2) THEN 
PI = 4.0 * ATAN(1.0) 


NAV = 3.0 
T=0 
H = 1.0 


XCENT = 4.0 
YCENT = 3.0 
SL = 200.0 
SB = 20.0 
IFLAG1(1) 
IFLAG2(1) 
GO TO 20 
ENDIF 
ELSE IF(IDRW.EQ.0) THEN 
GO TO 999 
END IF 
GO TO 130 
999 STOP 
END 


il i 
O 


© 
COARARARERAERA A AAA RRA RARA RRA AA RA RAR RARA RRA RARA RARA RARA RARA RA ARA AR 


SUBROUTINE AUTO1 (CSE2,U2) 


© THIS SUBROUTINE CALCULATES THE BEARING TO LAUNCH THE CHAFF AND 
THE 


C FIRST LEG OF THE NEW COURSES TO MANEUVER THE SHIP 
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COMMON/TIME/ITTG, TTG 
COMMON/WD/GW, VW 
COMMON/ABCD/BRGM(60) ,R(60) ,SC, VS (60) , CSE(60) , BRGMT (60), 
*BRGTM(60) , RMT (60) 

COMMON/MISS/RG, THETA, AG, PI 

COMMON/CHF/BRGC, ISHOOT 


COMMON/REL/XR (60), YR(60),XRD(60), YRD(60),XRP(60), YRP(60), XRPD(60), 


*YRPD(60) 
COMMON/TURN/TR, CSEOFF 
REAL CSEOPT(4),TOPT(4),CSE2,CSEDIF(4),CDIF1,ACC,VSMAX 
REALXRP1,YRP1,XRP2,YRP2,XRPD1,YRPD1,XRPD2,YRPD2,XRPDD, YRPDD 
INTEGER TTG1,T 
REAL CDIST 
POSSIBLE CHAFF LAUNCH BEARINGS 
BRGC1 = BRGM(1) + 90 
BRGC2 = BRGM(1) - 90 
IF (BRGC1.GT.360) THEN 
BRGC1 = BRGC1 - 360 
ELSEIF (BRGC1.LT.0) THEN 
BRGC1 = BRGC1 + 360 
ENDIF 
IF (BRGC2.GT.360) THEN 
BRGC2 = BRGC2 - 360 
ELSEIF (BRGC2.LT.0) THEN 
BRGC2 = BRGC2 + 360 
ENDIF 
WRITE (*,*) 'BRGCl=",BRGC1 
WRITE (*,*) 'BRGC2:',BRGC2 
GC = GW + 180 
VW1 = VW/(SC*1.8) 
COURSE OPTIONS FOR THE FIRST LEG 
CSEOPT(1) = BRGM(1) 
CSEOPT(2) = BRGM(1) + 180 
COURSE OPTIONS FOR THE SECOND LEG 
CSEOPT(3) = BRGM(1) + 90 
CSEOPT(4) = BRGM(1) - 90 
DO 10 I=1,4 
IF (CSEOPT(I).GT.360) THEN 
CSEOPT(I) = CSEOPT(I) - 360 
ELSEIF (CSEOPT(I).LT.0) THEN 
CSEOPT(I) = CSEOPT(I) + 360 
ENDIF 
CONTINUE 
CALCULATION OF INITIAL VELOCITY OF RELATIVE VELOCITY 
TEMP1=(VS(1)*(SIN((CSE(1)-BRGM(1))*PI/180))) 
TEMP2= (VW1* (SIN((GW-BRGM(1))*PI/180))) 
XRPD(1) = TEMP1+TEMP2 
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TEMP3=(VS(1)*(COS((CSE(1)-BRGM(1))*PI/180))) 
TEMP4=(VW1*(COS ( (GW-BRGM(1))*PI/180))) 
YRPD(1) = TEMP3+TEMP4 

C COURSE DIFFERENCE BETWEEN INITIAL COURSE AND THE COURSE OPTIONS 


DO 20 I = 1,4 
CSEDIF(I) = ABS(CSE(1) - CSEOPT(I)) 
IF (CSEDIF(I).GT.180) THEN 
CSEDIF(I) = 360 - CSEDIF(I) 
ENDIF 
20 CONTINUE 
C TIME TAKEN TO TURN TO COURSE OPTIONS 3 & 4 
b — 
TOPT(3) = CSEDIF(3)/TR 
TOPT(4) = CSEDIF(4)/TR 
INITIAL VALUES OF VARIABLES INVOLVED 
VSMAX = 25/(1.8*SC) 
TTG1 = IFIX(TTG) 
ACC = VSMAX * TR * PI/180 
T=1 
CRAD = 100./SC 
CDIST = (R(1) * TAN(THETA*PI/180)) - CRAD 
C  PRE-COMPUTATION TO CALCULATE CHAFF LAUNCH BEARING 


CH CH 


XRP1 = CDIST 
YRP1 = 0.0 
XRP2 = -XRP1 
YRP2 = 0.0 
XRPD1 = XRPD(1) 
YRPD1 = YRPD(1) 
XRPD2 = XRPD(1) 
YRPD2 = YRPD(1) 
XRPDD = 0.0 
YRPDD = 0.0 


C IF COURSE IS ALTERED TOWARDS COURSE OPTION 3(TO THE RIGHT) 


DO 30 I = 1,TTG1 
IF (I.LT.INT(TOPT(3))) THEN 


XRPDD = ACC * COS(ABS(CSE(1) - BRGM(1)) * PI/180) 
YRPDD = -ACC * SIN(ABS(CSE(1) - BRGM(1)) * PI/180) 
ELSE 
XRPDD = 0.0 
YRPDD = 0.0 
ENDIF 
XRP1 = XRP1 + XRPD1*T + 0.5 * XRPDD*T**2 
YRP1 = YRP1 + YRPD1*T + 0.5 * YRPDD*T**2 


XRPD1 = XRPD1 + XRPDD 
= YRPD1 + YRPDD 
30 CONTINUE 
© IF COURSE IS ALTERED TOWARDS COURSE OPTION 4 (TO THE LEFT) 


DO 40 I=1,TTG1 
IF (I.LT.INT(TOPT(4))) THEN 


XRPDD = -ACC * COS(ABS(CSE(1) - BRGM(1)) * PI/180) 
YRPDD = ACC * SIN(ABS(CSE(1) - BRGM(1)) * PI/180) 
ELSE 
XRPDD = 0.0 
YRPDD = 0.0 
ENDIF 
XRP2 = XRP2 + XRPD2*T + 0.5 * XRPDD*T**2 


YRP2 = YRP2 + YRPD2*T + 0.5 * YRPDD*T**2 
= XRPD2 + XRPDD 

YRPD2 = YRPD2 + YRPDD 

40 CONTINUE 

E CALCULTAION OF CHAFF LAUNCH BEARING 


IF (ABS(XRP1).GT.ABS(XRP2)) THEN 
BRGC = BRGC2 
ELSE 
BRGC = BRGC1 
ENDIF 
C CALCULATION OF THE FIRST LEG OF THE NEW COURSES 


IF (YRPD(1).GT.0) THEN 
CSE2 = CSEOPT(1) 
ELSE 
CSE2 = CSEOPT(2) 
ENDIF 
WRITE (*,*)'”CSE2=",CSE2 


E COURSE DIFFERENCE BETWEEN INITIAL COURSE AND THE FIRST LEG 

C 5 DT 
CDIF1 = ABS(CSE(1) - CSE2) 

E CALCULATION OF TURN DIRECTION INDICATOR FOR THE FIRST LEG 

C = — — e — = — — 


IF (CDIF1.GE.180) THEN 
IF (CSE2.GT.CSE(1)) THEN 
U2 = -1.0 
ELSE 
U2 = 1.0 
ENDIF 
ELSEIF (CDIF1.LT.180) THEN 
IF (CSE2.GT.CSE(1)) THEN 
U2 = 1.0 
ELSE 
U2 = -1.0 
ENDIF 
ENDIF 
RETURN 
END 
C kkkkkkkžkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkxk 
SUBROUTINE AUTO2 (CSE3,U3,CSE2) 
C THIS SUBROUTINE CALCULATES THE SECOND LEG OF THE NEW COURSES 
C TO MANEUVER THE SHIP 
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COMMON/MISS/RG,THETA, AG, PI 
COMMON/CHF/BRGC, ISHOOT 


COMMON/REL/XR(60) , YR(60) , XRD(60) , YRD(60) , XRP(60) , YRP(60) , XRDP(60), 
*YRDP (60) 
COMMON/ABCD/BRGM(60) ,R(60) ,SC,VS(60) ,CSE(60) BRGMT(60), 
*BRGTM(60) ,RMT(60) 
REAL CSEOPT (4) , CSE3, BRGM 


CSEOPT(1) = BRGM(1) 

CSEOPT(2) = BRGM(1) + 180 
CSEOPT(3) = BRGM(1) + 90 
CSEOPT (4) = BRGM(1) - 90 


DO 10 I=1,4 
IF (CSEOPT(I).GT.360) THEN 
CSEOPT(I) = CSEOPT(I) - 360 
ELSEIF (CSEOPT(I).LT.0) THEN 
CSEOPT(I) = CSEOPT(I) + 360 
ENDIF 
10 CONTINUE 
BRGC1 = BRGM(1) + 90 
BRGC2 = BRGM(1) - 90 
IF (BRGC1.GT.360) THEN 
BRGC1 = BRGC1 - 360 
ELSEIF (BRGC1.LT.0) THEN 
BRGC1 = BRGC1 + 360 
ENDIF 
IF (BRGC2.GT.360) THEN 
BRGC2 = BRGC2 - 360 
ELSEIF (BRGC2.LT.0) THEN 
BRGC2 = BRGC2 + 360 
ENDIF 
C CALCULATION OF THE SECOND LEG OF THE NEW COURSES 


IF (BRGC.EQ.BRGC1) THEN 


CSE3 = CSEOPT(4) 
ELSE 

CSE3 = CSEOPT(3) 
ENDIF 


WRITE (*,*)'CSE3=",CSE3 
CALCULATION OF THE TURN DIRECTION INDICATOR FOR 
THE SECOND LEG OF THE NEW COURSES 
IF (CSE3.EQ.CSEOPT(3)) THEN 
IF (CSE2.EQ.CSEOPT(1)) THEN 
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U3 = 1.0 

ELSEIF (CSE2.EQ.CSEOPT(2)) THEN 
U3 = -1.0 

ENDIF 


ELSEIF (CSE3.EQ.CSEOPT(4)) THEN 
IF (CSE2.EQ.CSEOPT(1)) THEN 
U3 = -1.0 
ELSEIF (CSE2.EQ.CSEOPT(2)) THEN 


99 


U3 = 1.0 
ENDIF 
ENDIF 
RETURN 


END 
Chkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


E 


SUBROUTINE PATH 
C THIS SUBROUTINE OF PLOT-88 DISPLAYS THE PATHS FOLLOWED 
C By THE SHIP, MISSILE AND CHAFF 
COMMON/PRN/IOPORT, MODEL 
COMMON/TIME/ITTG, TTG 
COMMON/TRACK/ XM(60) , YM(60) ,XS(60) , YS(60) ,XC(60) , YC(60) 
COMMON /MINN/RMIN 
REAL XCP (55) ,YCP(55) 
CHARACTER*20 LS,IM,LC,YD 
CHARACTER*25 LRMIN 
YD = 'VARDS' 
LRMIN="MINIMUM RANGE TO SHIP“ 
C SHIFT OF CHAFF POSITION ARRAY TO PLOT CHAFF AFTER IT BLOOMS 


DO 10 I = 1, (ITTG-5) 


XCP(I) = XC(I+5) 
YCP(I) = YC(I+5) 
10 CONTINUE 
C TO PREVENT UNDEFINED VALUES IF ITTG-5 < 55 
e a> a => > emm em em em em, em vm em gem em em em em em em em em em em em em e em «> > > > ov ge em em em em e e em em emm 
DO 20 I=(ITTG-4),55 
XCP(I) = XCP(I-1) 
YCP(I) = YCP(I-1) 
20 CONTINUE 
LS 2 's' 
LM = 'M' 
LC = 'C' 


K = 2*ITTG/3 

CALL PLOTS(0,IOPORT,MODEL) 

CALL CURVE (XM,YM,-ITTG,0.1) 

CALL SYMBOL(XM(K),YM(K),0.20,LM,0.0,1) 

CALL CURVE (XS,YS,-ITTG,0.1) 

CALL SYMBOL(XS(1),YS(1),0.20,LS,0.0,1) 

CALL CURVE (XCP, YCP,-ITTG,0.1) 

CALL SYMBOL(XC(6),YC(6),0.20,LC,0.0,1) 

CALL SYMBOL(0.5,0.2,0.25,LRMIN,0.0,21) 

CALL NUMBER(5.5,0.2,0.25,RMIN,0.0,1) 

CALL SYMBOL(7.0,0.2,0.25,YD,0.0,5) 

CALL PLOT(0,0,999) 

RETURN 

END 
Ckkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

SUBROUTINE SCREN1 
C THIS PLOT-88 SUBROUTINE DISPLAYS THE INITIAL CONDITIONS 
C OF THE SHIP AND THE MISSILE 
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COMMON/PRN/IOPORT, MODEL 
COMMON/CHF/BRGC, ISHOOT 
CHARACTER*25 CL 
CL = “AUTO CHAFF LAUNCH BEARING' 
CALL PLOTS(0,IOPORT,MODEL) 
CALL SHIP(1) 
CALL MISSLE(1) 
CALL RESCLL(1) 
CALL WIND 
CALL TABLE1(1) 
IF (ISHOOT.EQ.1) THEN 
CALL SYMBOL(1.0,5.7,0.16,CL,0.0,25) 
CALL NUMBER(6.0,5.7,0.16,BRGC,0.0,-1) 
END IF 
CALL PLOT(0,0,999) 
RETURN 
END 
CAŘRŘRŘRRÁRÁRÁ AAA RA AK K 
SUBROUTINE SCREN2 
C THIS PLOT-88 SUBROUTINE DISPLAYS THE INSTANTANEOUS POSITIONS 
C OF THE SHIP, CHAFF AND THE MISSILE EVERY 5 SECONDS 
COMMON/PRN/IOPORT, MODEL 
COMMON/TIME/ITTG, TTG 
COMMON/FG/IFLAG1 (60) , IFLAG2 (60) 
L = ITTG/5 
DO 31 J=1,L 
K=Jx*5 
IF (IFLAG2(N).EQ.1) THEN 
K = J+1 
END IF 
CALL PLOTS (0, IOPORT, MODEL) 
CALL SHIP(K) 
CALL MISSLE(K) 
IF (IFLAG1(K).EQ.0) THEN 
CALL RESCLL(K) 
END IF 
CALL WIND 
C CHAFF NOT DISPLAYED FOR THE FIRST 5 SECONDS TO SHOW BLOOM TIME 
DELAY 


IF (J.GT.1) THEN 
CALL CHAFF (K) 
ENDIF 
CALL TABLE1 (K) 
CALL PLOT(0,0,999) 
31 CONTINUE 
RETURN 
END 
Ckkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkáěkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


SUBROUTINE SHIP(J) 
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C THIS SUBROUTINE CALCULATES THE COORDINATES OF THE POLYGON 
C REPRESENTING THE SHIP 
COMMON/TRACK/ XM(60) , YM(60) ,XS(60) , ¥S(60) ,XC(60) , YC(60) 
COMMON/MISS/RG, THETA, AG, PI 
COMMON/ABCD/BRGM(60) ,R(60) ,SC, VS (60) ,CSE(60) , BRGMT (60), 
*+ BRGTM (60) , RMT (60) 
COMMON/SHP/SL, SB 
REAL XSHIP(5),YSHIP(5) 
C ROTATES THE SHIP DEPENDENT UPON THE COURSE 


XSHIP(1)=XS(J) + (SL/2.)*SIN(CSE(J) *PI/180.)/SC 
XSHIP(2)=XS(J) + (SL/4.)*SIN(CSE(J)*PI/180.)/SC + 
*(SB/2.)*SIN((CSE(J)+90.)*PI/180.)/SC 
XSHIP(3)=XS(J) - (SL/2.)*SIN(CSE(J)*PI/180.)/SC + 
*(SB/2.)*SIN((CSE(J)+90.)*PI/180.)/SC 
XSHIP(4)=XS(J) - (SL/2.)*SIN(CSE(J)*PI/180.)/SC + 
* (SB/2.)*SIN((CSE(J)-90.)*PI/180.)/SC 
XSHIP(5)=XS(J) + (SL/4.)*SIN(CSE(J)*PI/180.)/SC + 
* (SB/2.)*SIN((CSE(J)-90.)*PI/180.)/SC 
YSHIP(1)=YS(J) + (SL/2.)*COS(CSE(J)*PI/180.)/SC 
YSHIP(2)=YS(J) + (SL/4.)*COS(CSE(J)*PI/180.)/SC + 
* (SB/2.) *COS((CSE(J)+90.)*PI/180.)/SC 
YSHIP(3)=YS(J) - (SL/2.)*COS(CSE(J)*PI/180.)/SC + 
*(SB/2.)*COS((CSE(J)+90.)*PI/180.)/SC 
YSHIP(4)=YS(J) - (SL/2.)*COS(CSE(J)*PI/180.)/SC + 
* (SB/2.)*COS((CSE(J)-90.)*PI/180.)/SC 
YSHIP(5)=YS(J) + (SL/4.)*COS(CSE(J)*PI/180.)/SC + 
* (SB/2.)*COS((CSE(J)-90.)*PI/180.)/SC 
C THIS PLOT-88 SUBROUTINE FILLS THE SHIP POLYGON 


CALL STFILL(1) 
CALL FILL (xship,yship,5) 
RETURN 
END 
COA de de e e e e A A AR RAR RA A ARA A RA A A A A A A A e A A A e A A e e e e A A e e A A RA A KR KR 
SUBROUTINE RESCLL(J) 
C THIS SUBROUTINE CALCULATES THE INSTANTANEOUS COORDINATES OF 
C THE RADAR RESOLUTION CELL 
COMMON/MISS/RG,THETA,AG,PI 
COMMON/ABCD/BRGM(60) ,R(60) ,SC, VS (60) , CSE (60) , BRGMT (60), 
*BRGTM(60) , RMT (60) 
COMMON/EQ/ XET(60) , YET(60) 
REAL XCELL(4) , YCELL(4) 
C SET UP THE RESOLUTION CELL OF THE MISSILE 


RMBD = BRGTM(J) + 180. 
AG = RMT(J) * TAN(THETA*PI/180.) 
ALPHA = 90. + RMBD 


XW1 = XET(J) + AG*SIN(ALPHA*PI/180.) 
YW1 = YET(J) + AG*COS (ALPHA*PI/180.) 
XW2 = XET(J) - AG*SIN(ALPHA*PI/180.) 
YW2 = YET(J) - AG*COS(ALPHA*PI/180.) 
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XCELL(1) = XW1 + RG*SIN(RMBD*PI/180.) 
VCELL(1) = YW1 + RG*COS(RMBD*PI/180.) 
XCELL(2) = XW1 - RG*SIN(RMBD*PI/180.) 
YCELL(2) = YW1 - RG*COS(RMBD*PI/180.) 
XCELL(3) = XW2 - RG*SIN(RMBD*PI/180.) 
YCELL(3) = YW2 - RG*COS(RMBD*PI/180.) 
XCELL(4) = XW2 + RG*SIN(RMBD*PI/180.) 
YCELL(4) = YW2 + RG*COS (RMBD*PI/180.) 


C DRAW THE RESOLUTION CELL 
CALL STFILL(3) 
CALL FILL (XCELL, YCELL,4) 
RETURN | 
END 
Ckkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 
SUBROUTINE WIND 
C THIS PLOT-88 SUBROUTINE DISPLAYS THE WIND DIRECTION LINE 
COMMON/MISS/RG,THETA,AG,PI 
COMMON/WD/GW, VW 
CHARACTER*20 LW 


a 


LW ="W" 
XCENT = 4.0 
YCENT = 3.0 

© DRAW WIND DIRECTION LINE 
TWBRG = GW + 180. 
XWDP1 = XCENT - 2.8*SIN(TWBRG*PI/180.) 
YWDP1 = YCENT - 2.8*COS(TWBRG*PI/180.) 
XWDP2 = XCENT - 2.0*SIN(TWBRG*PI/180.) 
YWDP2 = YCENT - 2.0*COS(TWBRG*PI/180.) 


CALL PLOT(XWDP1,YWDP1,3) 
CALL PLOTD(XWDP2,YWDP2,2) 
CALL SYMBOL(XWDP1,YWDP1,0.20,LW,0.0,1) 
RETURN 
END 
CARRER RE RAR RS 
SUBROUTINE CHAFF(J) 
C THIS SUBROUTINE CALCULATES THE COORDINATES THE COORDINATES 
C OF THE HEXAGON REPRESENTING THE CHAFF 
COMMON/TRACK/ XM(60),YM(60),XS(60),YS(60),XC(60), YC(60) 
COMMON/ABCD/BRGM(60),R(60),SC,VS(60),CSE(60), BRGMT(60), 
*BRGTM (60), RMT(60) 
COMMON/CFCL/PARTC, CRAD, RGC, AGC 
COMMON /MISS/RG, THETA, AG, PI 
REAL XCF(6) , YCF(6) 
C © DRAW CHAFF CLOUD 


XCF(1)=XC(J) + CRAD 

YCF(1)=YC(J) 

XCF(2)=XC(J) + CRAD * COS(60.*PI/180.) 
YCF(2)=YC(J) + CRAD * SIN(60.*PI/180.) 
XCF(3)=XC(J) - CRAD * COS(60.*PI/180.) 
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YCF(3)=YC(J) + CRAD * SIN(60.*PI/180.) 
XCF (4)=XC(J) - CRAD 

YCF (4)=YC(J) 

XCF(5)=XC(J) - CRAD * COS(60.*PI/180.) 
YCF(5)=YC(J) - CRAD * SIN(60.*PI/180.) 
XCF(6)=XC(J) + CRAD * COS(60.*PI/180.) 
YCF(6)=YC(J) - CRAD * SIN(60.*PI/180.) 
CALL STFILL(2) 

CALL FILL (XCF,YCF,6) 

RETURN 

END 


Ckikkkikkikikkkkkkkikkkkikkkikkkkikkikkkikkikikikkikkikkikikikikikkikkikkikikikikkkkikikkkkik 


C 
FIRE 


SUBROUTINE MISSLE(J) 


THIS PLOT-88 SUBROUTINE DISPLAYS THE MISSILE BEARING /LINE OF 


COMMON/EQ/ XET(60) , YET(60) 
COMMON/MISS/RG, THETA, AG, PI 
COMMON/TRACK/ XM(60), YM(60),XS(60),YS(60),XC(60), YC(60) 
COMMON/ABCD/BRGM(60),R(60),SC,VS(60),CSE(60), BRGMT(60), 


*BRGTM (60), RMT(60) 


COMMON/FG/IFLAG1(60),IFLAG2(60) 
CHARACTER*20 LM 

LM = 'M' 

XM1 = XET(J) + 2.5 * SIN(BRGTM(J)*PI/180.) 
VMI 2 VET(J) + 2.5 * COS(BRGTM(J)*PI/180.) 


DRAW THE MISSILE BEARING LINE 


IF (IFLAG1(J).EQ.0) THEN 
IF (RMT(J).GT.2.5) THEN 
CALL PLOT(XET(J),YET(J) ,3) 
CALL PLOT (XM1,YM1,2) 
CALL SYMBOL(XM1,YM1,0.20,LM,0.0,1) 
ELSE 
CALL PLOT(XM(J),YM(J),3) 
CALL PLOT(XET(J), YET(J),2) 
CALL SYMBOL(XM(J),YM(J),0.20,1M,0.0,1) 
END IF 
ELSE 
CALL SYMBOL(XM(J),YM(J),0.20,1M,0.0,1) 
END IF 
RETURN 
END 


Chkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkxkxk 


SUBROUTINE TABLE1(J) 


C THIS PLOT-88 SUBROUTINE DISPLAYS THE VALUES OF IMPORTANT 
VARIABLES 

E BOTH FOR THE INITIAL PICTURE AND FOR THE RUN OF PLAY EVERY 5 
SECS 


COMMON/TME/T 
COMMON/ABCD/BRGM(60),R(60),SC,VS(60),CSE(60),BRGMT(60), 


*BRGTM(60) ,RMT(60) 


COMMON/WD/GW, VW 
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COMMON/SHCL/AGS , RGS, BRGMS (60) , SAA, PARTS, SAAN (60) 
COMMON/RC/RCSS ,RCSSP (60) , RCSC, RCSCP (60) 


C H A R A C T E R 
CS,SPD, WND,MLL,RGE,SCL,TR,AA,SRCS,CRCS,SP,BR,SCR, TM 
TM = våt 
BR = ‘’BRG’ 
SP = ‘SHIP’ 
CS = 'CSE” 
SPD = ‘SPD’ 
MLL = ”MISSILE' 
RGE = ‘RG’ 
WND = “WIND” 
SCR = 'SCREEN' 
SCL = ‘SCALE’ 
TR = ‘TRUE’ 
SRCS = ‘SHIP RCS’ 
CRCS = ‘CHAFF RCS’ 
C CONVERT BACK TO THE UNSCALED RANGE & SPEED 
C CE 
ZJ = J 
RR = R(J) * SC 
SHSPD = 1.8 * SC * VS(J) 
C DRAW TABLE AT RIGHT OF THE SCREEN 
C -- 
CALL PLOT (6.9,6.0,3) 
CALL PLOT (6.9,0.0,2) 
call symbol(7.0,5.8,0.18,SP,0.0,4) 
SAS mbol(7.0,5.5,0.16,CS,0.0,3) 
call number(7.6,5.5,0.16, CSE(J), 0.0,-1) 
call symbol(7.0,5:3,0.16,SPD,0.0,3) 
call number (7.6,5.3,0.16, SHSPD, 00 1) 
call symbol(7.0,4.8,0.16,WND,0.0, 4) 
call symbol(7.0,4.5,0.16,BR,0.0, 3) 
Gall number(7.6,4.5,0.16.6W,0.0,-1) 
call symbol(7.0,4.0,0.16, SPD, OF 0, 3) 
call number (7.6,4.0,0.16,VW,0.0, -1) 
call symbol(7.0,3.5,0.16,MLL,0. 0, 7) 
call symbol (7.0,3.0,0.16,BR,0.0, 7) 
call number (7.6,3.0,0.16, BRGM (J) , 0.0,-1) 
call symbol(7.0,2.5,0.16,RGE,0.0,2) 
call number(7.4,2.5,0.16,RR, 0. O, L1) 
call symbol(7.0,1.7,0.16,SCR,0. 0, 6) 
call svmbol(7.2,1.4,0.16,SCL,0.0,5) 
carl number(7.3,1.1,0.16,SC,0.0, 1) 
call symbol(7.1,0.8,0. 16,TM,0.0,3) 
call number(7.4,0.8,0.16,ZJ,0.0,-1) 
© SHIP RCS 
C æ æ > > > mv ov ov 
call symbol(4.0,0.1,0.16,SRCS,0.0,9) 
call number(6.0,0.1,0.16,RCSSP(J),0.0,1) 
C CHAFF RCS 
C ————— 


LOS 


call symbol(4.0,0.3,0.16,CRCS,0.0,9) 
call number (6.0,0.3,0.16,RCSCP(J),0.0,1) 
RETURN 
END 


Ckkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


SUBROUTINE SHPCLL(J) 
C THIS SUBROUTINE CALCULATES THE POTION OF THE SHIP RCS IN THE 
CELL 
COMMON/MISS/RG, THETA, AG, PI 
COMMON/ABCD/BRGM (60) ,R(60) ,SC,VS(60) , CSE(60) , BRGMT (60), 
*BRGTM (60) , RMT (60) 
COMMON/SHCL/AGS, RGS, BRGMS (60) , SAA, PARTS, SAAN (60) 
COMMON/SHP/SL, SB 


IF (RGS.GT. (RG-SL/(2.0*SC))) THEN 
IF (BRGMS(J).LE.0.0) THEN 
BR = 90.0 - (BRGMS(J) * 180./PI) 
ELSE 
BR = 450.0 - (BRGMS(J) * 180./PT) 
END IF 
CALL ASPECT (J,BR, SAAR) 
S2R = ABS(RGS-RG) 
IF (SAAR.EQ.90.0) THEN 
S1R = 0.0 
ELSE 
S1R = S2R / ABS(COS(SAAR*PI/180.)) 
END IF 
END IF 
IF (AGS.GT.(AG-SL/(2.0*SC))) THEN 
IF (BRGMS(J).GT.BRGMT(J)) THEN 
IF (BRGMS(J).LE.0.0) THEN 
BA = 180.0 - (BRGMS(J) * 180./PI) 
ELSE 
BA = 540.0 - (BRGMS(J) * 180./PI) 
END IF 
ELSE 
IF (BRGMS(J).LE.0.0) THEN 
BA = 0.0 - (BRGMS(J) * 180./PI) 
ELSE 
BA = 360.0 - (BRGMS(J) * 180./PI) 
END IF 
END IF 
CALL ASPECT (J,BA,SAAA) 
S2A = ABS(AGS-AG) 
IF (SAAA.EQ.90.0) THEN 
S1A = 0.0 
ELSE 
S1A = S2A / ABS(COS(SAAA*PI/180.)) 
END IF 
END IF 
IF((RGS.GT. (RG-SL/(SC*2.0))).AND. (AGS.GT. (AG-SL/(SC*2.0)))) 
THEN 
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IF ((RGS.LT.RG) .AND. (AGS.LT.AG)) THEN 
IF (S1A.LT.S1R) THEN 


Sl = SIA 

SAA = 90 - SAAA 
ELSE 

S1 = SIR 

SAA = SAAR 
END IF 


ELSE IF ((RGS.GT.RG).AND.(AGS.GT.AG)) THEN 
IF (S1A.LT.S1R) THEN 


S1 = SIR 
SAA = SAAR 
ELSE 
S1 = S1A 
SAA = 90. - SAAA 
END IF 
ELSE IF ((RGS.GT.RG).AND.(AGS.LE.AG)) THEN 
S1 = SIR 
SAA = SAAR 
ELSE 
S1 = S1A 
SAA = 90.0 - SAAA 
END IF 


ELSEIF((RGS.GT. (RG-SL/(SC*2.0))).AND. (AGS.LE. (AG-SL/(SC*2.0)))) 


*THEN 

Sl — SIR 

SAA = SAAR 
ELSE 

S1 = S1A 

SAA = 90. - SAAA 
END IF 


SAREA = SL * SB 
SAAN(J) = SAA 
S1SC = S1*SC 
IF ((RGS.GT.RG).OR.(AGS.GT.AG)) THEN 
PARTS = (SAREA - .5*SAREA - SB*S1SC) /SAREA 
IF (S1SC.GT.(SL/2.0)) THEN 
PARTS = 0.0 
END IF 
ELSE IF (((RGS.EQ.RG).AND. (AGS.LE.AG)).OR. 
* ((RGS.LT.RG).AND.(AGS.EO.AG))) THEN 
PARTS = 0.5*SAREA 
ELSE 
PARTS = (0.5*SAREA + SB*S1SC) /SAREA 
IF (S1SC.GT.(SL/2.0)) THEN 
PARTS = 1.0 
END IF 
END IF 
RETURN 
END 
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Chkkhkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 
SUBROUTINE ASPECT (J, B, SAA) 


C THIS SUBROUTINE CALCULATES THE ANGLE BETWEEN THE SHIP’S HEADING 
C AND THE GATE BOUNDARY IT EXITS. IT IS USED IN SUBROUTINE SHPCLL 
COMMON/ABCD/BRGM(60),R(60),SC,VS(60),CSE(60), BRGMT(60), 
*BRGTM(60) , RMT (60) 
C DETERMINE SHIP ASPECT ANGLE (SAA) 
C wem gn em em em wm em em em em em em em e em em em ëm mt Oe Oe ew ewww wwe = 
C QUADRANT II 
C = «=> av an av an am am am am am 
IF ((90.0.LE.B).AND.(B.LT.180.0)) THEN 
IF ((0.0.LE.CSE(J)).AND.(CSE(J).LT.90.0)) THEN 
IF (B.GE.CSE(J)) THEN 
SAA = B - CSE(J) 
ELSE 
SAA = CSE(J) + 180.0 - B 
END IF 
END IF 
IF ((90.0.LE.CSE(J)).AND. (CSE(J).LT.180.0)) THEN 
IF (B.GE.CSE(J)) THEN 
SAA = B - CSE(J) 
ELSE 
SAA = CSE(J) - B 
END IF 
END IF 
IF ((180.0.LE.CSE(J)).AND. (CSE(J).LT.270.0)) THEN 
IF ((CSE(J)-B).LE.90.0) THEN 
SAA = CSE(J) - B 
ELSE 
SAA = B - (CSE(J) - 180.0) 
END IF 
END IF 
IF ((270.0.LE.CSE(J)).AND.(CSE(J).LT.360.0)) THEN 
IF ((CSE(J)-B).LE.180.0) THEN 
SAA = B - (CSE(J) - 180.0) 
ELSE 
SAA = (CSE(J) - 180.0) - B 
END IF 
END IF 
C QUADRANT I 
(© ase ee eee ae — GE 


ELSE IF ((180.0.LE.B).AND.(B.LT.270.0)) THEN 
IF ((0.0.LE.CSE(J)).AND.(CSE(J).LT.90.0)) THEN 
IF (B.GE.(CSE(J)+180)) THEN 
SAA = B - CSE(J) 
ELSE 
SAA = CSE(J) + 180.0 - B 
END IF 
END IF 
IF ((90.0.LE.CSE(J)).AND.(CSE(J).LT.180.0)) THEN 
IF ((B-CSE(J)).LE.90.0) THEN 
SAA = B - CSE(J) 
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ELSE 
SAA = CSE(J) + 180.0 - B 
END IF 
END IF 
IF ((180.0.LE.CSE(J)).AND.(CSE(J).LT.270.0)) THEN 
IF (CSE(J).GE.B) THEN 
SAA = CSE(J) - B 
ELSE 
SAA = B - CSE(J) 
END IF 
END IF 
IF ((270.0.LE.CSE(J)).AND.(CSE(J).LT.360.0)) THEN 
IF ((CSE(J)-B).LE.90.0) THEN 
SAA = CSE(J) - B 


ELSE 
SAA = B - (CSE(J) - 180.0) 
END IF 
END IF 
C QUADRANT IV 
C SSS SS LO DIDI 
ELSE IF ((270.0.LE.B).AND.(B.LT.360.0)) THEN 
IF ((0.0.LE.CSE(J)).AND.(CSE(J).LT.90.0)) THEN 
IF ((B-(CSE(J)+180)).LE.90.0) THEN 
SAA = B - (CSE(J)+180.) 
ELSE 
SAA = CSE(J) - B + 360.0 
END IF 
END IF 
IF ((90.0.LE.CSE(J)).AND. (CSE(J).LT.180.0)) THEN 
IF (B.GE. (CSE(J)+180.0)) THEN 
SAA = B - (CSE(J)+180.0) 
ELSE 
SAA = CSE(J) + 180.0 - B 
END IF 
END IF 
IF ((180.0.LE.CSE(J)).AND.(CSE(J).LT.270.0)) THEN 
IF ((CSE(J)+180.0-B).GE.90.0) THEN 
SAA = B - CSE(J) 
ELSE 
SAA = CSE(J) + 180.0 - B 
END IF 
END IF 
IF ((270.0.LE.CSE(J)).AND.(CSE(J).LT.360.0)) THEN 
IF (B.GE.CSE(J)) THEN 
SAA = B - CSE(J) 
ELSE 
SAA = CSE(J) - B 
END IF 
END IF 
C  QUADRANT III 
C ID E ae a 


ELSE IF ((360.0.LE.B).AND.(B.LT.450.0)) THEN 
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IF ((0.0.LE.CSE(J)).AND. (CSE(J).LT.90.0)) THEN 
IF ((B-360.0).LE.CSE(J)) THEN 
SAA = CSE(J) - (B - 360.0) 
ELSE 
SAA = (B - 360.0) - CSE(J) 
END IF 
END IF 
IF ((90.0.LE.CSE(J)).AND. (CSE(J).LT.180.0)) THEN 
IF ((CSE(J)-(B-360.0)).LE.90.0) THEN 
SAA = CSE(J) - (B-360.0) 
ELSE 
SAA = (B - 360.0) - (CSE(J) - 180.0) 
END IF 
END IF 
IF ((180.0.LE.CSE(J)) .AND. (CSE(J) .LT.270.0)) THEN 
IF ((CSE(J)-180.0).GE.(B-360.0)) THEN 
SAA = (CSE(J) - 180.) - (B - 360.0) 
ELSE 
SAA = (B - 360.) - (CSE(J) - 180.0) 
END IF 
END IF 
IF ((270.0.LE.CSE(J)).AND. (CSE(J) .LT.360.0)) THEN 
IF ((B-CSE(J)).LE.90.0) THEN 
SAA = B - CSE(J) 
ELSE 
SAA = CSE(J) - (B - 180.0) 
END IF 
END IF 
ELSE 
END IF 
RETURN 
END 
C L A A RARA ARA KAKA 
SUBROUTINE CHFCLL 
C THIS SUBROUTINE CALCULATES THE CHAFF RCS IN THE CELL 
COMMON/MISS/RG,THETA, AG, PI 
COMMON/CFCL/PARTC, CRAD, RGC, AGC 
COMMON/RC/RCSS,, RCSSP (60) ,RCSC, RCSCP (60) 
CAREA — PI * CRAD**2.0 
RIR — ABS(RGC — RG) 
RIA — ABS(AGC — AG) 
IF ((R1R.LT.CRAD).AND.(R1R.GT.0.0)) THEN 
A2R = ACOS(R1R/CRAD) 
R2R = CRAD*SIN(A2R) 
SECTR = A2R*CRAD**2.0 
TRIR = R1R*R2R 
END IF 
IF ((R1A.LT.CRAD).AND.(R1A.GT.0.0)) THEN 
A2A = ACOS(R1A/CRAD) 
R2A = CRAD*SIN(A2A) 
SECTA = A2A*CRAD**2.0 
TRIA = R1A*R2A 
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END IF 
IF (((R1R.LT.CRAD).AND. (R1R.GT.0.0)).AND. 
* ((R1A.LT.CRAD) . AND. (R1A.GT.0.0))) THEN 


A2A = ACOS(R1A/CRAD) 
A2R = ACOS(R1R/CRAD) 
R2A = CRAD*SIN(A2A) 
R2R = CRAD*SIN(A2R) 


SECTA = A2A*CRAD**2.0 
SECTR = AZR*CRAD**2.0 


TRIA = RIA*R2A 
TRIR = RIR*R2R 
END IF 


IF ((R1R.EQ.0.0).OR. (RIA.EQ.0.0)) THEN 
PARTC = 0.5 
IF ((R1R.EQ0.0.0).AND. (R1A.EQ.0.0)) THEN 
PARTC = 0.25 
END IF 
ELSE IF ((R1R.LT.CRAD).AND. (R1A.GT.CRAD)) THEN 
IF (RGC.LT.RG) THEN 
PARTC = (CAREA-SECTR+TRIR) /CAREA 
ELSE 
PARTC = (SECTR-TRIR)/CAREA 
END IF 
ELSE IF ((R1R.GT.CRAD).AND.(R1A.LT.CRAD)) THEN 
IF (AGC.LT.AG) THEN 
PARTC = (CAREA-SECTA+TRIA) /CAREA 
ELSE 
PARTC = (SECTA-TRIA) /CAREA 
END IF 
ELSE IF (((R1R.LT.CRAD) .AND. (R1R.GT.0.0)).AND. 
*((R1A.LT.CRAD) .AND. (R1A.GT.0.0))) THEN 
IF ((RGC.LT.RG).AND. (AGC.LT.AG)) THEN 
IF (R1R.LT.R1A) THEN 
PARTC = (CAREA-SECTR+TRIR) /CAREA 
ELSE 
PARTC 
END IF 
ELSE IF ((RGC.LT.RG).AND. (AGC.GT.AG)) THEN 
PARTC = (SECTA-TRIA)/CAREA 
ELSE IF ((RGC.GT.RG).AND. (AGC.LT.AG)) THEN 
PARTC = (SECTR-TRIR) /CAREA 
ELSE 
IF (R1R.LT.R1A) THEN 
PARTC = 0.5*(SECTA-TRIA) /CAREA 
ELSE 
PARTC 
END IF 
END IF 
ELSE 
PARTC = 0.0 
END IF 
IF ((RGC.GT. (RG+CRAD)).OR. (AGC.GT. (AG+CRAD))) THEN 


(CAREA-SECTA+TRIA) /CAREA 


0.5* (SECTR-TRIR) /CAREA 


1il 


PARTC = 0.0 
END IF 
RETURN 
END 


Ckhkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


SUBROUTINE SCREN3 (IFLAG3) 


© THIS SUBROUTINE ENABLES ALTERATION OF THE INPUT VARIABLES 


100 


COMMON/ABCD/BRGM(60),R(60),SC,VS(60),CSE(60),BRGMT(60), 


*BRGTM (60), RMT(60) 


COMMON /WD/GW, VW 
COMMON/MISS/RG,THETA,AG, PI 
COMMON/RC/RCSS ,RCSSP (60) ,RCSC,RCSCP (60) 


INTEGER I 

IFLAG3=0 

WRITE (*,*) “(1) MISSILE BEARING’ 

WRITE(*,*) '(2) MISSILE RANGE’ 

WRITE(*,*) ’(3) SHIP’’S HEAD’ 

WRITE(*,*) "(äi SCREEN SCALE’ 

WRITE(*,*) "(Si MISSILE RADAR BEAMWIDTH’ 

WRITE(*,*) "(ei RANGE GATE’ 

WRITE(*,*) ”(7) MAXIMUM RADAR CROSS SECTION OF SHIP’ 
WRITE(*,*) "(8) MINIMUM RADAR CROSS SECTION OF SHIP“ 
WRITE(*,*) “(9) RADAR CROSS SECTION OF CHAFF’ 
WRITE(*,*) ’(10) TRUE WIND DIRECTION’ 

WRITE(*,*) “(11) TRUE WIND SPEED’ 

WRITE(*,*) “(12) ALL PARAMETERS 

WRITE(*,*) “(13) RETURN TO MAIN PROGRAM’ 


WRITE(*,*) INPUT SELECTION’ 
READ (*,*) I 
IF (I.EQ.1) THEN 
WRITE(*,*) ‘INPUT MISSILE BEARING’ 
READ (*,*) BRGM(1) 
IFLAG3 = 2 
GO TO 100 
ENDIF 
IF (I.EQ.2) THEN 
WRITE (*,*) ‘INPUT INITIAL MISSILE RANGE’ 
READ (*,*) RANGE 
IFLAG3 = 2 
GO TO 100 
ENDIF 
IF (I.EQ.3) THEN 
WRITE(*,*) ‘INPUT SHIP’’S HEAD” 
READ (*,*) CSE(1) 
IFLAG3 = 2 
GO TO 100 
ENDIF 
IF (I.E0.4) THEN 
WRITE(*,*) “INPUT SCREEN SCALE' 
READ (*,*) SC 
IFLAG3 = 2 
GO TO 100 
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ENDIF 
IF (I.EQ.5) THEN 
WRITE (*,*) ‘INPUT MISSILE RADAR BEAMWIDTH’ 
READ (*,*) THETA 
IFLAG3 = 2 
GO TO 100 
ENDIF 
IF (I.EQ.6) THEN 
WRITE (*,*) ‘INPUT RANGE GATE’ 
READ (*,*)RGATE 
RG=RG/SC 
IFLAG3 = 2 
GO TO 100 
ENDIF 
IF (I.EQ.7) THEN 
WRITE(*,*) ‘INPUT MAXIMUM RADAR CROSS SECTION OF SHIP’ 
READ (*,*) RCSSMAX 
IFLAG3 = 2 
GO TO 100 
ENDIF 
IF (I.EQ.8) THEN 
WRITE(*,*) ‘INPUT MINIMUM RADAR CROSS SECTION OF SHIP’ 
READ (*,*) RCSSMIN 
IFLAG3 = 2 
GO TO 100 
ENDIF 
IF (I.EQ.9) THEN 
WRITE(*,*) “INPUT RADAR CROSS SECTIONOF CHAFFf' 
READ (*,*) RCSSC 
IFLAG3 = 2 
GO TO 100 
ENDIF 
IF (I.EO.10) THEN 
WRITE(*,*) ‘INPUT TRUE WIND DIRECTION’ 
READ (*,*) GW 
IFLAG3 = 2 
GO TO 100 
ENDIF 
IF (I.EQ.11) THEN 
WRITE (*,*) ‘INPUT TRUE WIND SPEED’ 
READ (*,*) VW 
IFLAG3 = 2 
GO TO 100 
ENDIF 
IF (I.EQ.12) IFLAG3=1 
IF ((I.LT.1).OR.(I.GT.13)) THEN 
WRITE(*,*) “INPUT AGAIN’ 
GO TO 100 
ENDIF 
RETURN 
END 
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